Fluctuations and phase transitions in Larkin-Ovchinnikov liquid crystal states of 

population-imbalanced resonant Fermi gas 
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Motivated by a realization of imbalanced Feshbach-resonant atomic Fermi gases, we formulate 
a low-energy theory of the Fulde-Ferrell and the Larkin-Ovchinnikov (LO) states and use it to 
analyze fluctuations, stability, and phase transitions in these enigmatic finite momentum-paired 
superfluids. Focusing on the unidirectional LO pair-density wave state, that spontaneously breaks 
the continuous rotational and translational symmetries, we show that it is characterized by two 
Goldstone modes, corresponding to a superfluid phase and a smectic phonon. Because of the liquid- 
crystalline "softness" of the latter, at finite temperature the 3d state is characterized by a vanishing 
LO order parameter, quasi-Bragg peaks in the structure and momentum distribution functions, 
and a "charge"-4, paired Cooper-pairs, off-diagonal-long-range order, with a superfluid-stiffness 
anisotropy that diverges near a transition into a nonsuperfluid state. In addition to conventional 
integer vortices and dislocations the LO superfluid smectic exhibits composite half-integer vortex- 
dislocation defects. A proliferation of defects leads to a rich variety of descendant states, such 
as the "charge" -4 superfluid and Fermi-liquid nematics and topologically ordered nonsuperfluid 
states, that generically intervene between the LO state and the conventional superfluid and the 
polarized Fermi-liquid at low and high imbalance, respectively. The fermionic sector of the LO 
gapless superconductor is also quite unique, exhibiting a Fermi surface of Bogoliubov quasiparticles 
associated with the Andreev band of states, localized on the array of the LO domain-walls. 



I. INTRODUCTION 
A. Background 

1. Imbalanced resonant atomic gases 

Experimental progress in trapping, cooling and co- 
herently manipulating Feshbach-resonant atomic gases 
opened an unprecedented opportunity to study degen- 
erate strongly interacting quantum many-body systems 
in a broad range of previously unexplored regimes 
0). These include paired fermionic superfluids (SF) 

(with s-wave SF now readily realized iT Hlll . and p- 

wave SF under extensive current study [3|, I12h16I 1). the 
associated Bardeen-Cooper-SchriefFer (BCS) to Bose- 
Einstein condensation (BEC) crossover!^ [l7H25|. B ose- 
Fermi mixtures |26i] , bosonic molecular superfluids I27H3C 
and many other states and regimes [Slj under both equi- 
librium and nonequilibrium conditions [2^. |32| . Issj . In 
addition to the tunability of the Feshbach-resonant in- 
teraction strength (and its effective sign), temperature, 
and many types of external perturbations, a species num- 
ber imbalance in e.g., a two atomic hyperfine-states mix- 
ture turned out to be an extremely fruitful experimental 
knobfU-H^. 

A nonzero species imbalance frustrates conventional 
BCS pairing of a two-species Fer nii g as fssl - lilj and the 
associated BCS-BEC crossover |42h46| | . driving quantum 
phase transitions out of a paired superfluid to a variety 
of interesting pos sible ground states and thermodynamic 
phases [fl I46l45l| . This rekindled considerable theoreti- 
cal activity in the context of species-imbalanced resonant 
Fermi gases [1, [52| - [72| . The identification of the num- 
ber species imbalance with the magnetization of an elec- 



tronic system, and the chemical potential difference with 
an effective Zeeman energy, connects these recent atomic 
gases studies with a large body of research on solid state 
electronic superconductors under a Zeeman field [tI, [tsI - 
pTZj . as well as extensivelystudied realizations in nuclear 
and particle phvsics (78i - i81| . The obvious advantage of 
the current atomic system is the aforementioned tun- 
ability, disorder-free "samples", and absence of the or- 
bital part of the magnetic field, that always accompa- 
nies a solid-state charged superconductor in a magnetic 
field [821]. In these neutral paired superfluids the orbital 
field effects can be independently controlled by a rotation 
of the atomic cloud 18311. 



As illustrated in Fig|TJ among many interesting 
features, such as the gapless imbalanced superfluid 
fSf A.f) [4i - l47j , ubiquitous phase separationjUiO,!!! 47 1, 
tricritical point [sl,!?^, etc., observed experimentally 34- 
113] and studied extensively theoretically [1, 0], the 
interaction-imbalance BEC-BCS phase diagram is also 
predicted 0, H H HI IzO, llll] to exhibit the eiiigmatic 
Fulde-Ferrell-Larkin-Ovchinnikov state fFFLO)[76l. \7T\. 
First predicted in the context of solid-state superconduc- 
tors over 45 years ago|76i,i77,], the FFLO states has so far 
eluded a deflnitive observation, though some promising 
solid statefi^] and quasi-ld atomicHH candidate systems 
have recently been realized. 

At its most generic level the FFLO state[86] is a 
fermionic superfluid, paired at a flnite center of mass 
momentum. Generically such a state spontaneously 
"breaks" gauge and translational symmetry, i.e., it is a 
periodically-paired superfluid (superconductor), akin to 
a supersolid[87-_9Q], and thus can appropriately be called 
a pair-density wave (PDW)[9ll l9^ . Microscopically, it 
is driven by Fermi surface mismatch (tgI. [77| due to an 
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FIG. 1: A mean-field zero-temperature phase diagram from 
Refsi, m of an imbalanced Fermi gas, as a function of tlie 
inverse scattering length and normalized species imbalance 
P = [N^-Ni)/N = AN/N, showing the magnetized (imbal- 
anced) superfluid (SFm), the FFLO state (approximated as 
the simplest FF state, confined to a narrow red sliver bounded 
by Pfflo and Pc2) and the imbalanced normal Fermi liquid. 

imposed pairing species number (and/or massfosf) im- 
balance. As a compromise between the superfluid pair- 
ing and an imposed imbalance, at intermediate values 
of the latter, the superconducting order parameter con- 
denses at a set of finite center-of-mass momenta deter- 
mined by the details of the Fermi surface mismatch and 
interactions, thereby self-consistently leading to FFLO 
pairing between these imbalanced fermionic species. At 
sufficiently large imbalance, AiVc2 or equivalently at the 
upper-critical Zeeman field hc2 (corresponding to the 
chemical potential difference Afic2 of the two pairing 
species) no compromise is possible, and instead a transi- 
tion to the normal state takes place. 



2. Mean-field energetic stabtUty of FFLO states 

The original predictions by Fulde and Ferrell (FF), and 
by Larkin and Ovchinnikov (LO) were followed by ex- 
tensive studies of the FFLO states [78-8 though ex- 
clusively confined to their energetics in the BCS limit. 
Stimulated by the aforementioned potential realization 
in imbalanced resonant Fermi gases, the recent revival of 
the subject extended the analysis to the full range of the 
BCS-BEC crossoverdHlllSllO,!!!. As illustrated 
in Fig 111 the key observation is that, despite of strong 
interactions, within a single qo BdG treatment the con- 
ventional FFLO state (those originally proposed by FF 
and LO) remains quite fragile, confined to a narrow sliver 
of polarization in the BCS regime [3. |4^. 

However, motivated by earlier studies of the 
Bogohubov-de Gennes (BdG) eauation [95l - l97| . combined 
with finding of a negative domain-wall energy in an 
otherwise fully-paired singlet BCS superfluid in Zeeman 
field [13 , recent studies have quite convincingly ar- 
gued, that a more generic pair-density wave state (that 




FIG. 2: An illustration of a continuous commensurate- 
incommensurate (CI) transition at hci from a fully- 
gapped (balanced) paired-superfiuid to an imbalanced Larkin- 
Ovchinnikov superfluid. The excess majority atoms are local- 
ized on the domain walls in (zeros of) the LO order parameter, 
whose number ndw(h) is then proportional to the imbalance 
P{h) and grows continuously with the chemical potential dif- 
ference (Zeeman energy), h — hc\. 



includes a larger set of coUinear wavevectors) may be sig- 
nificantly more stable. Analogously to a strongly type- 
11 superconductor, that undergoes a continuous transi- 
tion in to a vortex state at a lower-critical orbital field 
i/c 1 [9^-[l0il (that is significantly below the thermody- 
namic field, He) in the current system, a Zeeman-field 
h (the chemical potential imbalance) can drive a nu- 
cleation of domain-walls in the superfluid order pa- 
rameter above the lower-critical Zeeman field hd that 
is below the bulk mean- field value he- Thereby, the 
excess of the majority fermionic atoms (polarization) 
in an imbalanced system can be continuously accom- 
modated by the sub-gap states localized on the self- 
consistently induced domain- walls, through a continuous 
commensurate- incomme nsur ate (CI) Pokrovsky-Talapov 
(PT) type of transition lo^l from a fully paired s-wave 
superfluid to a Larkin- O vchinnikov like periodic state 
of domain walls [l03l Il04| . This picture, illustrated in 
Figl2] resemb les t he soliton mechanism for doping of 
polvacetvlene jlOSj l . The ±A domain- wall description, 
that is explicitly realized in one- dimension (Id) through 
exact BdG [9^ ar id B e the ansatz (l06lll07| solutions and 
via bosonization |l08l Il09l |. is complementary, but not 
qualitatively distinct from the more familiar single co- 
sine form of the LO state [tJ. 
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The latter, "soft" pair-density wave is a more appro- 
priate description near the hc2 transition into the normal 
state, where the pair-order parameter is naturally small, 
and a Landau expansion in the leading harmonic, is 
expected to be valid. On the other hand, clearly, far be- 
low hc2 (e.g., near hd transition from the fully-paired 
uniform singlet BCS state) the periodic soliton state is 
a quantitatively better description as it more accurately 
captures the strong pairing at A, confining the "normal" 
regions that accommodate the imposed fermion imbal- 
ance to the narrow gapless domain-walls between ±A. 
Upon increasing h above hd the density of domain- walls 
grows, eventually overlapping at hc2, thereby interpo- 
lating between the two limiting descriptions of the LO 
state. This picture is explicitly realized in the exact Id 
solutions [95I lOd ^ 108 1 and emerges from the numerical 
BdG studies [91193" Furthermore, the lack of a single q 
nesting for d > 1, to which the fragility of the LO state 
is usually attributed is irrelevant when the LO order pa- 
rameter exhibits a broad spectrum of q (set by l/^o), as 
in its soliton form above hd- 

The LO state can be equivalently thought of as a peri- 
odically ordered micro-phase separation between the nor- 
mal and paired states, that naturally replaces the macro- 
phase separation [13, |8l| ubiquitously found in the BCS- 
BEC detuning-imbalance phase diagram 

SB US- 

It is clear that the standard (even multi-wavevector) 
Landau treatments valid near hc2 and analytical BdG 
analysis of a single FF plane-wave state (exhibiting no 
amplitude nodes) Q fail to capture above quantitatively 
important ingredients. These treatments are therefore 
not necessarily quantitatively trustworthy in their predic- 
tion of the energetic range of stability (location of hd) of 
the LO state (they are however reliable for the prediction 
of hc2), and in our view need to be reexamined. 



3. Fluctuation tn the FFLO states 



fermionic atoms confined to a fl uctuating periodic array 
of two-dimensional domain walls jllll Ill2| . 

In the isotropic realization (e.g., in cold atoms in 
an isotropic trap) of interest to us, the FFLO states 
spontaneously break a continuous rotational (in addi- 
tion the translational) symmetry akin to smectic liq- 
uid crystals, in contrast to their solid state density-wave 
analogs. C onse quently, as was originally anticipated by 
Shimaharajnoi, and was demonstrated in our recent 
workfTT^I, to be explored in greater detail below, their 
Goldstone modes are qualitatively "softer" , and there- 
fore exhibit far stronger fluctuations. These can either 
completely destabilize the (otherwise energetically sta- 
ble) FFLO state, or can qualitatively modify its mean- 
fleld form and properties. This general picture therefore 
reveals that the complexity of the FFLO state beyond 
its mean-fleld cartoon requires the understanding of sub- 
tle interplay of superfluidity, liquid-crystallization, and 
anisotropic Fermi surface physics. 

Wit h th e above motivation in mind, in a recent 
Letter [11 25 above questions were formulated and care- 
fully explored. Namely, as illustrated in a schematic 
phase diagram in Figl21 supported by the afore mentioned 
studiesfil Hi] and the exact Id solutionsfH [lol [Tol], 
we assumed that a striped (unidirectional) FFLO state 
(that we refer to as LO), characterized by a collinear 
set of wavevectors is energetically (microscopically) sta- 
ble over an experimentally accessible portion of the 
detuning-imbalance phase diagram. We then formulated 
the model for the low-energy Goldstone modes fluctua- 
tions and fermionic excitations in the LO state and used 
it to study the stability of the state to quantum and 
thermal fluctuations, as well as explored the fluctuation- 
driven phenomenology, topological defects, novel quan- 
tum liquid crystal phases, and the associated phase 
transitions [11 2j . In the present manuscript we present 
a significantly more extensive description of these find- 
ings and the details of the associated calculations. 



The microscopic question of the energetic stability of 
FFLO states is certainly an extremely important one and 
has dominated most of the research on the subject to 
date. However, given the extensive 45-year history of 
the topic, it is astounding that the equally basic comple- 
mentary question of the nature of Goldstone modes de- 
scription and their fluctuat ions within the FFLO states 
received so little attention, jllOl . Illl| until our study of 
the problem, reported in a recent Letter [112|. From 
the general symmetry principles the FFLO states' low- 
energy phenomenology is expected to be signiflcantly 
richer than that of a homogeneous fully gapped supercon- 
ductor, whose low-energy phenomenological (Ginzburg- 
Landau and xy-model) desc riptio n long predated the mi- 
croscopic theory by BCS [9i-[l0l|. Namely, in addition to 
a local superfluid phase degree of freedom, the low-energy 
modes include the phonons (a single scalar one in the 
case of a uniaxial LO state) of the periodic superconduct- 
ing structure, as well as gapless polarized (single-species) 



4-. Relation to other systems and studies 

Although physicall y qu ite distinct, some of our mo- 
tivation and findings |ll2j | are closely related to stud- 
ies in solid state realizations, such as the putative 
FFLO states in heavy- fermion (CeCoIns) a nd organi c 
(k- (BEDT-TTF)2Cu(NCS)2 ) superconductors [ul [ill , 
the striped (and spiral) states in high-Tc superconduc- 
tors and nickelates and more generally strongly cor- 
related doped Mott in sulators [11 5- 12(11, spiral states 
in helimagnets [l2ll4l23l. par tially-filled high Landau 
level 2d electron gas 124| - |l26l |. charge-density waves in 
anisotropic metals 127| . and others. Among a number 
of physical properties special to resonant atomic Fermi 
gases, one key qualitatively distinguishing feature in our 
study is that the FF and LO states spontaneously breaks 
continuous translational and orientational symmetries of 
an isotropically trapped Fermi gas. As first emphasized 
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FIG. 3: A proposed AN/N vs. l/(fcFas) phase diagram for 
an imbalanced resonant Fermi gas, showing the more stable 
LO liquid crystal phases (discussed in the text and illustrated 
in detail Fig|4} replacing a portion of the phase-separated 
regime. 

in Refs llldlll2l this latter property is responsible for the 
qualitatively enhanced fluctuations, that lead to the pre- 
dicted universal power-law correlations (l5]),([71) persisting 
throughout the 3d quasi-long-range ordered FFLO phase, 
rather than just at a critical point. Thus, such FF and 
LO states (and its descendents) are qualitatively distinct 
from their mean- field form[76l, [77|. akin to a distinction 
between, e.g., a 2d xy model and a long-range ordered 
3d ferromagnet. In contrast, the aforementioned states 
in the solid state systems, by their very definition only 
break discrete point-group crystal symmetries, and there- 
fore exhibit quite tame Gaussian fluctuations inside the 
phase, with these periodic states only quantitatively dis- 
tinct from their mean-field form. 

However, as we will see, despite enhanced fiuctuations 
the LO state studied here does lead to a number of 
novel features recently discussed in the literature, such 
as the fractional vortex-dislocation defects, "charge" - 
4 superflui dity, Andreev-like mid-gap states, and many 

others tnilnl]. 



5. Outline 

The rest of the paper is organized as follows. We con- 
clude the Introduction with the summary of our main 
results. Then, in Sec|lT]we use a microscopic description 
of the imbalanced resonant atomic Fermi gas to derive in 
the weakly-interacting BCS limit the corresponding Lan- 
dau theory (along with microscopic expressions for its 
parameters) valid near hc2 transition to the normal state. 
In Sec lIIII we discuss symmetries of two generic class of 
FFLO states, construct corresponding order parameters 
and derive from the Landau theory the corresponding 
models for their Goldstone modes. As we argue, these 
models provide universal descriptions of low-energy fluc- 
tuations in the FF and LO class of states, beyond the 



limit of validity of their microscopic BCS-limit deriva- 
tion. In Sec lIVI we complement this derivation (valid 
near hc2) by an analysis of an array of fluctuating domain 
walls, valid near hd- In SeclVl^e use these low-energy 
universal models to analyze quantum and thermal fluctu- 
ations in the FF and LO class of states, paying particular 
attention to elastic nonlinearities that (as we show) are 
qualitatively important at a finite temperature. Then in 
Sec|VT] we classify and discuss energetics of topological 
defects in the LO state. In Sec lVIII we use this infor- 
mation to uncover a variety of exotic "daughter" liquid- 
crystal phases that emerge as a result of unbinding dif- 
ferent combinations of topological defects, and discussion 
corresponding phase transitions. In Sec I Villi we extend 
the model to include the coupling of the Goldstone modes 
to gapless fermionic atoms, and comment on their effects. 
In Sec lIXI we use local density approximation (LDA) to 
put our results in the context of trapped atomic gases and 
discuss experimental probes of our predictions. We con- 
clude in Secs lXIXIIXIII with discussions of experimental 
probes, open questions and a summary of our study. In 
appendices A and B we provide technical details for the 
microscopic derivation of the Ginzburg-Landau expan- 
sion near tic2 and finite-size scaling analysis of smectic 
fluctuations, respectively. 

B. Results 

Before turning to the derivation and analysis of 
the model, we summarize the kevOTedictions of our 
work, previously reported in Ref ill 21 . Because our 
predictions 128] are based on general symmetry princi- 
ples, supported by detailed microscopic weak-coupling 
calculations, they are generic and robust to variation in 
microscopic details. At a very general level, we demon- 
strate that (in contrast to the conventional uniform su- 
perfluid and FF states) a unidirectional (striped) LO su- 
perfluid exhibits two Goldstone modes, </> and 9. These 
correspond to two coupled smectics phonons, or equiva- 
lently the superfluid phase (p and a nonlinearly-coupled 
smectic phonon mode u = —O/qq, with qo the wavevector 
characterizing the LO state. 

Through robust symmetry arguments, complemented 
by an explicit microscopic derivation (valid in the BCS 
regime, near the LO to normal state transition at hc2), we 
show that the low-energy universal (classical) LO Hamil- 
tonian governing fluctuations of these Goldstone modes 
is given by 

nfl>' = |(vM2-K|(a|i^-i(v«)^)%|(v.0)^ 

(1) 

This Hamiltonian form, familiar from studies of conven- 
tional smectic liquid crystals, encodes the underlying ro- 
tational invariance through the vanishing of the (Vj^u)^ 
modulus and the specific form of the nonlinear elastic 
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terms in n |l29| -fT32!|. In above ||, _L refer to axes that are 
parallel and perpendicular to the LO ordering wavevector 
qo, respectively. In the weak-coupling BCS limit, we de- 
rive explicit expressions for the above Goldstone-modes 
moduli, K, B, pi {i =|j,_L), given by 



epqo 



(2a) 



B « ^2^1 ^ 3.3nA|^^ Hh/h^,), (2b) 



(2c) 



also given in Eas. ([5n|) . ([5T|) of the main text. Among 
these predictions, we find that the LO state is a highly 
anisotropic superfluid, with the ratio of superfluid stiff- 
nesses given by 

pi/pl = J (A,„/Abcs)' « I Hhc2/h) « 1, (3) 



vanishing on the approach to the upper-critical Zeeman 
field h — >■ h~2, that marks the mean- field transition to 
the normal state at which the LO order parameter Ag^ 
vanishes. Abcs is the zero-field BCS order parameter. 

We find the FF state to be even more exotic. In con- 
trast to other homogeneous superfiuids (described by an 
xy- model), its single superfluid Goldstone mode cj) is de- 
scribed by above smectic Hamiltonian, with an identi- 
cally vanishing transverse superfluid stiffness, pj: pp = 0, 
a reflection of the rotational invariance of the sponta- 
neous current to an energy-equivalent ground state. Thus 
we show that a resonant imbalanced Fermi gas, confined 
to an isotropic trap gives a natural realization of a quan- 
tum (superfiuid) smectic liquid crystal. 

As a consequence of the spontaneous breaking of con- 
tinuous spatial symmetries (contrasting with its solid 
state realizations, where only discrete point-g roup (lat- 
tice ) symmetries are spontaneously broken) |113l 111 



TT8l4mfl23l - [T25t the Goldstone mode excitations in the 
FF and LO states are of qualitatively lower energy. As a 
result, the fluctuations in such superfiuid smectic states 
are qualitatively stronger. Specifically, we find that while 
they are stable to quantum fiuctuations, in 3d the LO 
and FF long-range orders are marginally unstable at any 
nonzero T, with the LO order parameter (with Dirichlet 
boundary conditions; for a more general case, see Sec 1X1 
of the main text) 



(Aio(r))fl, = (2A 



90 



e*'*('') cos (qo • r - 



e{v))), 



— cos qn • r 



0, 



(4) 



vanishing in the thermodynamic limit (a large cloud with 
atom number and cloud size i? — > oo), suppressed to 
zero by thermal phonon u = —O/qo fiuctuations, and is 
therefore strictly speaking homogeneous on long scales. 
The resulting superfluid state is thus an "algebraic topo- 
logical" phase with no long-ranged translational order. 



Namely, beyond mean-field theory it is instead character- 
ized by power-law order-parameter correlations, distin- 
guished from the spatially short-ranged disordered phase 
by confined topological defects (bound dislocations), not 
by a nonzero LO order parameter. It is therefore a 3d 
analog of the more familiar quasi-long-range ordered su- 
perfluid fllm and a 2d easy-plane ferromagnet. 

There are a number of interesting consequences of this 
finding. F or ex ample, as illustrated in Fig|2Tl we predict 
that in 3d |l33l | the static structure function S'(q) in the 
LO state exhibits universal anisotropic gttasi- Bragg peaks 
(akin to a conventional smectic liquid crvstal jl34 Il35j |). 
with n-th order peak given by 



^(9|hq± = 0) 



\q\\ - 2?igop"^"'''^ 



(5) 



rather than the true i5-function Bragg peaks of e.g., a 
long-range-ordered pair-density wave in a crystallin e en- 
vironment. In above the anomalous Caille exponent (l34l | 
is given by 



(6) 



We similarly find that the momentum distribution 
function of pairs displays a power-law form around the 
reciprocal lattice momenta set by of q^, 



1 



nqa 



2 — n^7] 



(7) 



as illustrated in Fig. (j20p . This power-law is a re- 
flection of a striking pair-condensate depletion to zero 
by the divergent finite T LO Goldstone-mod e fluctua- 
tions even in 3d, akin to the Landau- Peierls |l36l Il37j | 
behavio r of film s of a conventional superfluid and 2d 
cr ystalsflsMill- Such static correlations in the LO 
state can be computed asymptotically exactly, as was 
first done for a conventional smectic liquid crvstal jl3ll 
I132I Il34| . This fluctuation-driven 3d power- law phe- 
nomenology is a unique feature of a unidirectional 
[collinear wavevectors) FFLO state. It is not exhibited 
by crystalline FFLO phases with multiple non-collinear 
ordering wavevectors [tsI [79j . that, in contrast are charac- 
terized by the long-range positional order and a nonzero 
pair-condensate, that is stable to thermal fluctuations. 

As with treat ments of t he LO state, where long-range 
order is assumed |ll3Llll8j |. in this algebraic LO phase we 
also find an unusual topological excitation - a half vortex 
bound to a half dislocation - in addition to integer vor- 
tices and dislocations. These are illustrated in Figs l8l9[ 
andllOl 

In 2d, at nonzero T the LO state is even more strongly 
disordered, at intermediate scales characterized by uni- 
versal power-law phonon correlations and concomitant 
short-range positional order with Lorentzian structure 
function peaks, coi itrol led by a nontrivial exactly cal- 
culable fixed point ^13. Asymptotically, however, at 
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arbitrary low temperature the 2 d LO state is unsta- 
ble to proliferation of dislocations jl42| . The state that 
results from such dislocated superfluid smec tic is ei- 
ther a "c harge" -4 (paired Cooper pairs) Il4^ nematic 
superfluid jll2l Ill9j | or a nema tic (possibly "fractional- 
ized"b) Fermi liquid [llTl Il26j |. latter qualitatively the 
same as the deformed Fermi surface state [49l|. 

More generally, while analyzing defects-driven contin- 
uous transitions out of the LO state, we uncover a rich 
array of descendent states, that generically must inter- 
vene between the LO superfluid and a fully-paired con- 
ventional (isotropic and homogeneous) superfluid and a 
conventional Fermi liquid. If indeed the 3d LO state is 
energetically stable, as argued above, we expect these 
novel states to appear in the region collectively denoted 



"LO liquid crystals" of the detuning-polarization phase 
diagram of FiglSl They include a nonsuperfluid smectic 
{FL^g^, driven by an unbinding of integer 27r- vortices) , 
and a superfluid {SF^, driven by a proliferation of inte- 
ger a-dislocations) and a nonsuperfluid {FL^, driven by 
an unbinding of both vortices and dislocations) nemat- 
ics, and the corresponding isotropic states, when disclina- 
tions also condense. In addition, we predict a variety of 
topologically-ordered isotropic and nematic "fractional- 
ized" Fermi-hquid states {FL*j^, FL*j^, FL}, and others), 
that are distinguished from their more conventional fully- 
disordered forms by gapped (bound) half-integer defects. 
These phases are summarized by a schematic phase dia- 
gram illustrated in FiglH We now turn to the derivation 
of these predictions. 




SF, (BCS) 

FIG. 4: A schematic imbalance-chemical potential (Zeeman energy), h = jii^ — m vs detuning (interaction strength), —l/kpa 
phase diagram, illustrating the 3d LO smectic phase (SFsm) and its descendant (described in the text), driven by a proliferation 
of various combinations of topological defects. The inset shows the global imbalance-interaction BCS-BEC phase diagram, 
illustrating the location of these putative phases. 



II. BCS THEORY OF IMBALANCED 
RESONANT FERMI GAS 

A. Model 

We begin with the one-channel model of two-species, 
resonantly interacting Fermi gas, appropriate for the ex- 



perimentally relevant broad Feshbach resonance [^-|6[. In 
the grand-canonical ensemble, it is described by a Hamil- 
tonian 

k,(T 
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with the two atomic species (labeling the hyperfine states 
=tj i) open-channel fermions created by the anticom- 



muting operator c^^ 



kcr' ^k'o-'} 



k,k' 



(9) 



with the single-particle energy — fi,^fc^/2m, mass m, 
and system's volume V (henceforth set to 1). 

In above we introduced two distinct chemical poten- 
tials, /Ltcr = {iJ,^,fii) to impose numbers of two separately 
conserved atomic species, No- = (iV^, iV^), or equivalently 
the total fermion number = A'^ -I- and the atom 
species imbalance A = — . Equivalently, the two 
species chemical potentials fj,^ = ^ + h and = j.i — h 
are related to the total-number chemical potential /i and 
the species-imbalance chemical potential h, latter cor- 
responding to the pseudo Zeeman energy. The imbal- 
anced resonant Fermi gas thermodynamics as a func- 
tion of N, AN, T, as, i.e., the extension of the BEC-BCS 
crossover to a finite imbalance can be computed by a va- 
riety of theoretical techniques, including quantum Monte 
Carlo [12, mean-field theorWj,!!!!!!,!!!!!!!, the large- 
Nf (fermion flavor) |l44 , Il45j | and e-expansions [l46j | . 

The attractive interaction is parameterized by a short- 
range s-wave pseudopotential with a strength g < . 
Through a standard T-matrix scattering calculation jl47l 
Il48l |. that gives the two-atom s-wave scattering ampli- 
tude, fsik) — —j^Tk, the pseudo-potential parameter 
g can be related to th e exp erimentally determined, mag- 
netic field dependent jl49| | scattering length 



1 1 V- 1 



(10) 



where the ultraviolet-divergent second term is regularized 
by a microscopic momentum cutoff scale A ^ 1/d set by 
the range of the potential, d. This gives (with h = 1 
hereafter) 



, , /47r 2A 

asig) = h — 

\mg IT 

m g 



A-K 1 + glg^ 



(11a) 
(lib) 
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where g^ is the effective screened coupling and g^ = 
is the critical coupling g (set by the zero-point energy at 
scale d) at which a molecular bound state appears and 
the scattering length diverges. The above relation allows 
a definition of the model and therefore a reexpression 
of physical observables in terms of the experimentally 
defined (UV-cutoff independent) scattering length Ug- 

To treat the many-body problem, we utilize the 
standard mean-field analysis d, HI] (quantitatively valid 
deep in the BCS regime, kpag ^ 1, but expected to 
be qualitatively valid throughout) by first assuming an 
expectation value 



5(c;(r)ct(r)) = A(r), 



(12) 



corresponding to pair condensation at a superposition 
of finite momenta q, with the set of Aq and q to be 
self-consistently determined. With this mean-field as- 
sumption, H, in Eq. ([5]), reduces to the standard BCS 
mean-field form: 



H 



q ^ k 



(13) 



i:(A 

q.k 



q'^k+a^c_k^a^ 



'^'-k+ft'^k+a^^q 



that can equivalently be obtained using a Hubbard- 
Stratonovich transformation and a saddle-point approx- 
imation on the coherent state path-integral formulation 
of the problem. 

Although the resulting mean-field Hamiltonian is 
quadratic in the fermionic operators, its diagonalization 
for a generic A(r) (an arbitrary set of Fourier components 
Aq) is only possible through a numerical self-consistent 
solution of the Bogoliubov-de Gennes equations [oBl [96j. 

Analytical progress is however possible through two 
complementary approaches. One is to specialize to a 
single Fourier component, Aq FF state[7^, that is self- 
consistently determined through the q-dependent gap 
equation, (jl2p . equivalent to the ground-state energy 
minimization. This approach allows for a computa- 
tion of the ground-state energy that is fully nonlinear 
in Aq. However, as emphasized in the Introduction 
with regard to learning about a more generic amplitude- 
modulated FFLO state this simplifying specialization 
to a single Fourier component is only harmless near 
the normal-to-FFLO state transition at the mean-field 
upper-critical field /ic2, where Aq is small and physics is 
well-approximated by the lowest dominant harmonic. 

An alternative analytical approach is through the 
Ginzburg-Landau expansion, that allows a treatment of 
a general form of the superconducting order parameter, 
(|12l) . but is again limited to the vicinity of hc2, where 
A(r) is small, permitting Taylor expansion of the ground- 
state energy in terms of itf77|. 



Bogoliubov— de Gennes analysis of the 
Fulde-Ferrell state 



1. Ground state solution 

Specializing to the Fulde-Ferrell (single harmonic, Aq) 
state simplifies the mean-field Hamiltonian, (jSL allow- 
ing for a straightforward diagonalization |4l |44|. [Ta] of the 
the associated Bogoliubov-de Gennes equations. A stan- 
dard analysis 4] gives the (fermionic) Bogoliubov quasi- 
particle operators, dka 



akt 



UkC_ 



-J 



a 



(14a) 
(14b) 
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with the orthonormal coherence factors [99l-[l0H 



with 



1- 

V2\ 



(15a) 
(15b) 



ensuring the canonical anticommutation relations, 

{"^kcr' '^kcr'} ~ ^<y,<j'- In terms of these quasi-particles, 
the Hamiltonian reduces to a diagonal quadratic form 

Hff = J2 (£^k.e(£;k.)aL«k. - EkaQ{-Ek^)a^^al^ 



kcr 



+ i?£s(Aq), 

with the ground state energy given by 



(16) 



(e^+A^)l/^ 



(19a) 
(19b) 



Egs 



k 



(efc - Ek + Skte(-£;kt) + E^^ei-E^^)) 



(17) 



In above, we defined the excitation energy iJko- for spin 
state a 



E, 



kt 



E, 



ki 



Eu 



E. 



2m ' 

kg 

2m ' 



(18a) 
(18b) 



and /i = ^(/^^ + /^^), h = ^(/i^ — IJ>\,) the average and 
difference chemical potentials, respectively. The intro- 
duction of 0-functions in above expressions divides the 
momentum space into three regions, for h > Q given by: 
(1) k e ki such that E]^-^ > and Ei^i > 0, (2) k e k2 
such that i?kt < and (3) k S ks such that Ei^i < 0, 
(it is easy to see that the E^g- cannot both be negative) , 
and by construction ensures a positive-definite excitation 
spectrum jii'kcrl above the FF ground state 



iFFr. 



n 



n 



''k+fi 

kgk3 kek2 



n ^-k+Ht n ("k + ^'k4+a/_k+at)|0)' 



(20a) 
(20b) 



keki 



with |0) the atom vacuum. The first form demonstrates 
that the FF ground state can be formally thought of 
as the Fermi sea of negative energy Bogoliubov "quasi- 
particles" added to the BCS state \BCSq) (that is not 
the ground state) with the center of mass momentum q. 
The second form, above, shows that equivalently the FF 
ground state corresponds to a finite center of mass mo- 
mentum (q) BCS state with unpaired atoms from the 
two momenta sets k2 and ks, defined above. 



This analysis thus suggests a possibility of two distinct 
FF states, FFj^i and FFys2, that exhibit a single- and 
two-species Fermi surfaces, respectively, with their vol- 
ume difference set by the imposed species imbalance [621]. 



T^FF 

Egs 



2. Analysis near hc2 

The energetic stability of the state is controlled by 
(Aq), that deep in the BCS limit was first computed 
by Fulde and Ferrellj?^. Its full behavior throughout the 
BCS-BEC crossover generally requires a combination of 
numerical and analytical analysis. This led to the predic- 
tion of a narrow sliver of stability of the FF state confined 
to the BCS side of the crossover, -l/(fcFas) > 0.46jl|4ll. 
However, as discussed in the Introduction, given the ex- 
pected unfavorable nature of the FF state (lacking the 
amplitude nodes necessary to accommodate the imbal- 
anced atoms) this analysis is unlikely to shed light on 
the stability of other (e.g., LO type) states far below the 
hc2 transition to them. Nevertheless, because the transi- 
tion at the upper-critical field hc2 (at least in mean-field 
theory) is expected to be continuous and therefore Aq 
small, growing from 0, E^g (Aq) can be analyzed ana- 
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lytically by Taylor-expanding it in Aq. 

Taking advantage of the extensive analysis in[3, IS], 
near hc2 we obtain 



E, 



FF 
GS 



£9|Aq| 



where 



3n 



1 . Vpq'^ - 4/i^ 



-In 

2 



and 



'-^BCS 



3n 



in 



q,q,q,q 



ieF vlq-^ - 4/i2 ' 



(21) 
2h 

~2h\ ' 
(22) 

(23) 



that agree with the Larkin-Ovchinnikov result [77[ ob- 
tained via a direct Landau expansion in Aq. In above 



we used (jlOp to eliminate the interaction coupling g in 
favor of the scattering length as and then reexpressed 
the latter in terms of the BCS superfluid gap A bcs ac- 
cording to ~m/{'i'Kas) = A^(eF) In (Se^^/Ascs), with 
3d density of states N{e) = m^/"^ ^/e/ Near its 
minimum Sq can be approximated quadratically in q 



3n 

4e_F 



In 



2a- 



Vp 

1 



2Va^ 



intf 



^BCS 



in 


■ h ' 


= - — In 






hc2_ 



(24) 

(25a) 
(25b) 

c.- - 1 IQklql ■ ^^^^^ 

In above, a is a solution of equation arising from the 
minimization of the ground-state energy with respect to 
center of mass momentum grj^, |4^ 

a — -In -, (26a) 



with 
£0 = 

90 = 
J s= 



1 



1 , a 
- In — 

2 a 

1.200. 



1' 



(26b) 

Combining this value with the point at which £o(^) van- 
ishes, we obtain the N-FFLO upper-critical transition 
fieldfzilll 



^FFLO 



^BCS 



« 0.754Abcs 
At this transition point we have 
a ^Bcs 



(27a) 
(27b) 



9o = 



Abcs 0.58 



1.81- 



3n 



Co 



1 16ej;^ ql 



n 1 

0.61 J, at the hc2 transition 

eF qo 



(28a) 

at the hc2 transition, (28b) 
(28c) 
(28d) 



where the BCS coherence length is given by its standard 
expression = vp / {■k/S.bcs)- 

Near the hc2 transition the FF quartic vertex, (|23|) also 
simplifies giving 



3n 



1 



1 



q,q,q,q 



Aep - 1 4/i2 ' 
3n 0.57 



Aep 

3n 1 



Aep A2 



(29a) 
(29b) 

at the hc2 transition, (29c) 



BCS 



which agrees precisely with perturbative LO result j77|. 

As we will see shortly, in single plane-wave FF state 
the transverse superfluid stiffness for a supercurrent flow- 
ing transversely vanishes identically, enforced by a Ward 
identity. As mentioned in the Introduction this prop- 
erty is a special feature of a single momentum com- 
ponent FF state and is guaranteed by the underlying 
rotational invariance of the spontaneous current to an 
energy-equivalent ground state. Thus, to correctly (even 
qualitatively) capture a more general FFLO (e.g., the 
LO) state requires the inclusion of the multiple momen- 
tum components, as in and the analysis beyond the 
lowest-order LO treatment [t^ • This unfortunately can- 
not be calculated analytically through the Bogoliubov- 
de Gennes analysis, although it quite successfully can 
nume ricallyjH ^ and through the Ginzburg-Landau 
expansion to which we now turn. 



C. Ginzburg-Landau expansion near hc2 

The analytical treatment of the LO state near hc2 relies 
on the Ginzburg-Landau expansion in Aq, that is small 
near the (in m ean-fi eld) continuous hc2 normal-to-FFLO 
transition [t^, Il50l |. This expectation is explicitly sup- 
ported by the exact Id BdG solution[95| at high fields, 
where A (a;) is indeed well-approximated by a single har- 
monic, with an amplitude Aq that vanishes continuously 
near hc2- 

Based on these general arguments, near hc2 the 
Ginzburg-Landau expansion for the ground-state energy 
is expected to take a familiar form 



n 



2^^,|Aq| 



E 



V A* A A* A 

''qi ,q2 ,q3 ,q4 ^qi ^q2 ^q., ^q4 7 



qi,q2,q3,q4 



(30) 



where rotational invariance constrains Eq to be a function 
of the magnitude of q only, q4 = qi — q2 + Qs , and Aq 
is a Fourier transform of A(r). 

It is straightforward to see that translational and ro- 
tational invariances and the quadratic form (in Aq) of 
the first term in H guarantee that Eq is independent of 
the type of the FFLO state, and is therefore identical to 
that of the FF state, ([2^ 4]. This observation guarantees 
that all harmonics with magnitude qo become unstable at 
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the same imbalance field hc2, (j27bp . with the degeneracy 
only lifted by the quartic Aq term in H. 

Following the standard prescription and guided by 
the seminal (lowest order) analysis of Larkin and 
Ovchinnikov^TTj, we now derive the quartic vertex func- 
tion, V'qj,q2,q3,q4, appearing in Eq. (pO|) . It, together with 



Eq will then allow us to derive the key elastic moduli 
(compression modulus, bending rigidity, and superfluid 
stiffnesses) characterizing the Goldstone mode s of the 
striped FFLO states. As first found in Ref ill 2| . one key 
observation is that a derivation of a generic Ginzburg- 
Landau form and consequently of a generic Goldstone 
mode theory requires a calculation that is of higher 
order than that originally carried out by Larkin and 
Ovchinnikov[77j . In particular, as we will see below, the 
leading momentum dependence of the quartic coupling 
needs to be kept. It corresponds to an induced short- 
ranged current-current interaction, Vijjijj, with the su- 
percurrent given by the standard expression 



J. = -Re[-A*(r)ia,A(r) 

TO 



(31) 



To derive the generic Ginzburg-Landau form, we 
use the imaginary time fermionic coherent-state path- 
integral formulation of the BCS problem, by computing 
the partition function. 



where for notational short-hand we defined a space-time 
coordinate x = (r, r) and a single particle Hamiltonian 

operator, I<t = - ^ - Ma- 

We use the Cooper-pair Hubbard-Stratonovich field 
A(x) to decouple the (quartic) pairing interaction g and 
to integrate out the fermionic atoms: 



[dipadip^dAdAle 



[dA*dA] 



,-Srff[A*,A] 



(34a) 
(34b) 



where Sr ^ So + S^nt - I^P- 



5-0 



1, 



S^nt = / [Ar^i>l + iJii^^A*], 



(35a) 
(35b) 



and we defined the Ginzburg-Landau effective action 



S'eff[A*,A] = -In 



dipcrje' 



-Srli'l,ll>^.A',A] 



(36) 



Trace 



(32a) 
(32b) 



where H is the four- Fermi Hamiltonian, ?A(j(t, r) 
Grassmann (anticommuting) fields, and the imaginary 
time action is given by 



5^ ['(/'*, V'tr] = 



drd'^r 



(33) 
I 



where 



Taylor-expanding Sr in powers of Smt we compute ^off 
in powers of A. Relegating all technical details to Ap- 
pendix |A] and focusing on the time- independent quar- 
tic order (in A(r)) contribution, given by the connected 
fourth cumulant ^4 = J drH^ , we find 



i?4 = 



y(ri,r2,r3,r4)A*(ri)A(r2)A*(r3)A(r4), 

(37) 



y(ri,r2,r3,r4) = / ^G?(r2 - n, c^)G;(r2 - ra, -c^)G°(r4 - rg, w)G°(r4 - ri, -c^). 



V'(qi,q2,q3,q4) = 



'-■6, ^j^^y^A 

^|-p^GO(k,c.)G°(qi - k, -w)G?(k - qi + q2,c.)G;(q4 - k, -c.). 



(38a) 
(38b) 



with the noninteracting fermionic Green's function (in We focus on the FFLO type order parameter 
Fourier space) as usual given by 

A(r) = ^Aq„(r)e''i"-, (40) 
q„ 

Q /, ,*\ ,on \ with Fourier transform given by 

G^{uJn,q) = -(■0<T^o-)o, (39a) 

= , I . (39b) ^(q) = E^q"(q~q")- (41) 
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We note that to go beyond mean-field, in above we in- 
cluded an additional long-scale positional dependence in 
the Larkin-Ovchinnikov order parameters A^^^ (r) on top 

I 



of the short-scale mean-field periodic dependence at the 
LO wavevectors Qn encoded in the plane-wave factor. 
Substituting this form into H4, gives 



-q-.4.o 



(42) 



xV'(q„i +qi,qn2 +q2,qn3 +q3,qn4 +q4)A*_^^(qi)Aq„Jq2)A*^Jq3)Aq„^(q4) 

I 



with <\i = — q„; ^ q„. (due to A^^, (q^) expected to 
be sharply peaked around q^ = 0), allowing us to dis- 
entangle the q sums and integrals, and the correspond- 
ing 5-function above, without double-counting momen- 
tum states. 

To lowest order in q^, we ignore its dependence in 
V, equivalent to the LO treatment [xlj, with t^(q„i -I- 

qi,qn2 + q2,q«3 +q3,qn4 + q4) ~ "^(qni,q™2,qn3,q«4)7 

and thereby obtain the quartic vertex that determines 



J 



which set of reciprocal momenta q„; (satisfying the mo- 
mentum conservation) minimizes the interaction energy, 
thereby defining the structure of the FFLO ground state. 
As found by LOjT^, near hc2 it is the striped state of 
coUinear q„; 's that is energetically favored and is well ap- 
proximated by the lowest pair of harmonics, q„; = ±q. 
Focusing on this energetically preferred (near hc2) LO 
state reduces the general form of in Eq. (^5)) to 



1 



Hi = -j (2^)'^(5'*(qi-q2+q3-q4) 



^ y(qi +qi,qi +q2,qi +q3,qi +q4)A*^(qi)Aq,(q2)A*^(q3)Aqi(q4) 
qi=±q 



+4F(q + qi, q + q2, -q + qa, -q + q4)A;(qi)Aq(q2)Alq(q3)A_q(q4) 

I 



(43) 



Taylor-expansion of F(q„i -1-qi, q„2 +q2, qnj +q3, qn4 
q4) in qi then gives = Hf^ + with 



z;(°|(|Aq(r)r + |A_q(r)r 



4t;(°l|Aq(r)nA_q(r)|^ 



where 

,(0) 



y(q, q, q,q). 



(44) 



(45a) 



/(^G?(k,.)^G°(q-k,-.) = 



w+1 = V'(q,q, -q, -q), 

dojd'^k ~ 



J (^G°(k,.)^G;(q-k,-.)G;(-q-k, 



-w 



already computed by LO in their mean-field approxima- 
tion and in Ref. ^ from the expansion of the BdG ground 

state energy of the FF state discussed in Sec lII Bl 

(2) 

In computing the higher order HI term, that contains 
the essential current-current interaction Vijjijj discussed 



in the previous section, we will neglect contributions from 
the first set of terms in Ea. (|43|) (-I--I- and terms), be- 
cause they lead to (positive q current) and j^q (neg- 
ative q current) contributions whose coefficients are en- 
forced either by rotational invariance (which for trans- 
verse current each are 1/2 of the coefficient of the jqj-q 
term, required to keep the smectic phonon field u "soft" ) , 
or for the longitudinal pieces are small, i.e., higher order, 
subleading corrections to the nonzero quadratic (in A) 
terms computed in the previous section. 



To proceed we focus on the second term in Eq. (|43 
and Taylor expand it to second order in q.; 



(45b) F(q- 



qi,q + q2, 

„(2), 



q3, 



q4) 



(46) 



+ 



'(q) 



4m'^ 



taking advantage of momentum conservation qi — q2 + 
q3 — q4 = to generate four equivalent terms such that 
the expression is explicitly real. In above we defined the 
current-current coupling matrix 
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,(2) 



(47a) 
(47b) 



where 

91 



\VFqoJ 



are computed in Appendix [X] Converting this Fourier 
space expression into real space gives 



(48a) 



2.07 



N{eF)kj 



92 



N{eF)kl 



-5.75 



Q!2 



at the hc2 transition, (48b) 
■ (48c) 



N{eF)P 



vfQo 

F 



at the hc2 transition, (48d) 



(2) 



2 / 



Irc (-A**a.A,(r)) Irc (-Al,^9, A_q(r)) 



2 J «J JqJ-q' 



24'^(jq+j-q).(jq+j-q) 



(49a) 
(49b) 
(49c) 
(49d) 



where we reconstructed the j^jl and jL^jL^ transverse 
pieces (by rotational invariance) and defined the total 
supercurrent j = (jq + j_q)/2. 

It is important to note already at this stage that the 
positivity of gi > in Eq. (|48|) generates a positive trans- 
verse superfluid stiffness, and therefore a well-defined 
LO state. Furthermore, we note that (as we will see in 
Sec lIIip the fact that 172 and even gi + 172 are negative 
does not cause any stability difficulties as they lead to 
0[{Aq„/ AbcsY] corrections to the longitudinal super- 
fluid stiffness pi, a correction that is subdominant near 
hc2- However, the opposite signs do suggest that below 
/ic2, P't grows, while pi decreases with decreasing chemi- 
cal potential difference, h. It is thus conceivable that the 
anisotropy ratio p^/pi' (starting at just below hc2) may 
actually grow above 1, i.e., that the superfluid anisotropy 
may reverses, something that would have striking exper- 
imental consequences. 

Putting above derived ingredients together we finally 



obtain the sought-after Ginzburg-Landau Hamiltonian 

Hgl - J [iV^Ap - 2g2|VA|2] +r|Ap-f iAi|A|4+iA2j'+. . • , 

(50) 

where as just derived, deep in the BCS limit (large pos- 
itive detuning, fci^a^ ^ 1) and near the hc2 transition 
to the polarized normal state the model parameters are 
given by 



J R 


0.61n 

4 ' 
epqo 




(51a) 


qo - 


l.SlAscs 

hvF ' 


(51b) 


r R 


3n , 
i in 

AeF 


■ 9h ' 

4/lc2. 


(51c) 


hc2 « 


6 . 

i -Abcs, 


(51d) 


Ai « 


3n 

" "^^F Ales' 


(51e) 



J^.LJKfl 111 I U 

efAIcs^o 
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Based on the discussion following Eq. we approxi- 
mated the anisotropic tensor coupling v^^ , Eq. (jTZj) by an 

(2) 

isotropic one of strength set by the transverse part of v^j , 
namely, gi, with the difference a subleading correction 
near hc2- More generally (away from the weak-coupling 
BCS limit) these couplings can be taken as phenomeno- 
logical parameters to be determined experimentally. 

III. THEORY OF GOLDSTONE MODES IN 
STRIPED FFLO STATES 

A. Landau theory 

The general form of the Ginzburg-Landau model ([50)) 
derived in the previous section has a much broader 
range of applicability, even if our derivation and micro- 
scopic predictions for the associated coupling constants 
in Eq. (|5ip only apply in the BCS regime, near hc2- Inde- 
pendent of the microscopies, the key ingredient of Ugl 
is that it captures the imbalanced atomic Fermi system's 
energetic tendency to pair at a finite momentum, and 
thereby forms a pair-density wave characterized by a re- 
ciprocal lattice vector with magnitude qo and sponta- 
neously chosen orientation. More generally, we expect 
the quadratic (in A(r)) part of Hgl 

= A*eA, (52) 

to be characterized by a more generic differential kernel 
(specialized to a quadratic form, e = J(— — q^Y + f' — 
JqQ in (jSOp ). exhibiting a minimum at a finite momentum 
with a magnitude go and an arbitrary orientation. How- 
ever, it can be shown jl5l| that this generalization adds 
little new physics to the Lifshitz-like normal-to-FFLO 
(PDW) transition and to the emerging Goldstone mode 
theory, that is our main interest. We will there work 
directly with the simpler form in (jSOp . 

With the minimum in the dispersion e{q) located at a 
finite momentum magnitude qo, for r < (that in the 
mean-field BCS approximation is given by rc — JqQ ~ 
0.61n/eF, or equivalently for at /i < hc2 « j'^bcs) the 
gap in e{q) closes and Hgl in (f50| develops an instability 



where f = r — Jq^ + J{q^ — q^Y . The quadratic coupling 
f dictates that the most unstable momentum mode is 
q ~ qo, that condenses when f becomes negative, i.e., r 
falls below = Jqo, corresponding to h = hc2- Since 
the first and second terms are "rotationally invariant" in 



to a nonzero superconducting (pairing) order parameter 
^(i") = X)q ^q,iS*'*" a set of nonzero wavevectors 
q„ with a magnitude of the fundamental given by the 
dispersion minimum go- 

As with other crystallization problems, the nontrivial 
question of the choice of the set of momenta q„ 's is de- 
termined by the details of the interaction (terms higher 
than quadratic order in A) and will not be re-addressed 
here. Instead, motivated by the LO findings and by the 
more recent analyses [tJ l95l - [98| , we will focus on the uni- 
directional pair-density wave (Cooper-pair stripe) order, 
characterized by a collinear set of q„'s. That is, we will 
assume that such states are energetically stable, will de- 
velop their Goldstone modes low-energy description and 
analyze their stability to fluctuations. 

As we will see below, within this unidirectional pair- 
density wave class of FFLO states, the FF and LO states 
are representatives of two qualitatively distinct univer- 
sality classes and therefore must be treated separately. 
Although (as argued in the Introduction) it is the latter 
that is expected to be significantly more stable, for com- 
pleteness and potential of other microscopic realizations 
(where FF may be stable) we will treat both universality 
classes. 

The low-energy properties of the FF and LO states 
are described by a periodically spatially modulated order 
parameter A(r), that in its simplest form, quantitatively 
valid near hc2 is well captured with a single ±q pair of 
ordering momenta 

Afflo (r) = A+ (r)e^'i-"- + A_ (r)e-''^•^ (53) 

where A± (r) are two complex scalar order parameters, 
the dominant Fourier coefficients of A(r) 

A±(r) = A^(r)e''^±(''), (54) 

and amplitudes Aj. distinguishing between the FF and 
LO states. We first focus on the amplitudes of these 
two order parameters, for now ignoring the corresponding 
Goldstone modes 4>± . Taking A± (r) as spatially indepen- 
dent (to be justified a posteriori), we use the Ginzburg- 
Landau theory, (j50p to determine their magnitudes. To 
this end we find 



(55) 

I 

the A+ — A_ space, i.e., 0(4) invariant, it is the last 
term that breaks this symmetry down to the physical 
U{1) (g) C/(l) and thereby determines the relative size of 
these critical ±qo momenta order parameters, A±. 

Clearly, for Ai > X2:3-, ^Landau is minimized by the 



I 

ULandau = r (|A+p + |A_p) + \x, (|A+p + |A_p)^ + i(A2^ - (|A+P - |A_p)^ , 
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FF state, with only one of the two order parameters 
nonzero, with = and 

2 

In the opposite limit of Ai < X2^^, the last quartic 
term in ULandau instead selects the LO state, with the 
two order parameters equal and nonzero, 

^LO ^ ^LO ^ ^|f|/(3Ai), for Ai < (57a) 

= ^AlcsHh/hc2), (57b) 

^ K'^ih). (57c) 

As first found by LO[77| and discussed in the previous 
section, in the current system the microscopies dictates 
that it is the latter, LO state that is the more stable one. 

B. Symmetries and order parameters 

The fundamental FFLO order parameter, ([55)) . that 
characterizes the FF and LO states clearly distinguishes 
these two symmetry-distinct states. 

1. The Fulde-Ferrell state 

The FF state is characterized by a single (independent) 
nonzero complex order parameter, 

Aff{v) = A,„e''i°-+'^ (58) 

that is a plane- wave with the momentum qo and a single 
Goldstone mode 

= '^+, (59) 

corresponding to the local superconducting phase. The 
state carries a nonzero, uniform spontaneously-directed 
supercurrent 

jfF = l|A,j2(q„ + V(/)), (60) 

and thereby breaks the time-reversal and rotational sym- 
metry, chosen spontaneously along qo, as well as the 
global gauge symmetry, corresponding to the total atom 
conservation. Although the FF order parameter itself 
is not translationally invariant, under translation by an 
arbitrary vector a it transforms by a multiplication by 
a global phase e''^" '^. It is therefore invariant under 
a modified transformation of an arbitrary translation 
followed by a gauge transformation. Thus, in the FF 
state all gauge-invariant observables and therefore the 
state are translationally invariant. Namely, the FF state 
is a uniform orientationally-ordered (polar) superfluid. 



Under an infinitesimal rotation of the FF current axis 
Qo — ^ Qo + (^qo; its phase transforms as (/>—>■(/) -I- (5qo • r, 
costing zero energy. Thus the underlying rotational sym- 
metry of the FF state requires the corresponding Gold- 
stone mode Hamiltonian of the superconducting phase 
= to be invariant under such transformation, that 
for an infinitesimal rotation corresponds to a phase shift 
linear in r transverse to qo . A generic Goldstone- mode 
Hamiltonian, that is well-known to satisfy these proper- 
ties is that of a smecticfnlllli], to harmonic order given 

by 

= ^^(Vi0)^ + ip!(qo-V0)2. (61) 

Its key qualitative feature is the strict (symmetry- 
enforced) vanishing of the (Vj_0)^ stiffness, with _L desig- 
nating axes transverse to the spontaneous current (fy = 
z = qo) axis. Thus, despite its uniform density (i.e., it is 
an orientationally ordered superfluid^ rather than a den- 
sity wave) the FF state is characterized by a smectic-like 
Goldstone-mode Hamiltonian, that is qualitatively dis- 
tinct from that of a conventional uniform and isotropic 
superfluid, described by an xy-model, 

H.y = ^P.(V0)^ (62) 

Namely, as evident from Eq. (IFTI) the FF state is in- 
finitely anisotropic, characterized by an identically van- 
ishing transverse superfluid stiffness 

pZ = 0, (63) 

a reflection of its underlying rotational symmetry that is 
spontaneously broken by the ground-state supercurrent 
'}FF- Its longitudinal superfluid stiffness Ps^\\ is nonzero, 
measuring the energetic cost of a deviation from the 
ground-state current magnitude oi j%p — ^|Aq„pgo- 

In the next section (after discussing the LO state), we 
will support these symmetry-based arguments by an ex- 
plicit derivation from the generic Ginzburg-Landau the- 
ory dnni) for the FF state. 



2. The Larkin-Ovchinnikov state 

As illustrated through the Landau analysis, the LO 
state is described by a nonzero (standing- wave like) pair- 
density wave order parameter. That is, taking the mag- 
nitudes of the two order parameters, A-|_ = A_ = Aq^ 
in ([5^ to be the same (as dictated by the last quartic 

term for Ai < ^2-;^) the LO order parameter reduces to 
a physically appealing form 

Aio(r) = 2A,„e^5W++'A-)cos[qo-r + i(0+ -(/._)], (64a) 
= 2A,„e**cos[qo-r + 0], (64b) 
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that is a product of a superfluid and a unidirectional den- 
sity wave order parameters, respectively characterized by 
two Goldstone modes 



(65a) 
(65b) 



In a qualitative contrast to the FF state and to a conven- 
tional superfluid or a superconductor, LO state's density 
and other physical quantities are periodic along qg, ex- 
hibiting periodic uniaxial stripe order. The position of 
the associated pair-density wave is characterized by the 
smectic phonon 



u = -9/qo, 



(66) 



giving physical interpretation to the second Goldstone 
mode 9, as the phase of the pair-density wave. 

We also note that unlike a conventional smectic [l29j | 
(e.g., in liquid crystal materials, where one instead is 
dealing with a real mass density p(r) not a pair con- 
densate wavefunction) , here, because A(r) is complex, 
the phases of A+are independent (though interacting) 
Goldstone modes p^53 |. 

The mean-field LO order parameter, A^o thus si- 
multaneously exhibits the ODLRO (superfluid) and the 
smectic (unidirectional density wave) orders. It thus 
spontaneously breaks the rotational, translational, and 
global gauge symmetries, and is therefore realizes a form 
of a paired supersolid. However, it is distinguished from 
a conventional purely bosonic supersolid [STJooj, where 
homogeneous superfluid order and periodic density wave 
coexist, by the vanishing of the ("charge"-2 two-atom) 
zero momentum (q — 0) superfluid component in the LO 
condensate 91]. 

The supercurrent in the LO state is given by 



JLO 



2IA P 

' * '-V(/)[l + cos(2q-r + 26')], (67a) 



m 



(67b) 



where in the last form we neglected its periodic contri- 
bution. As expected, in contrast to the FF state, (1501) . 
the supercurrent vanishes in the LO ground state where 
V^^ 0. 

Similarly to the FF state, the underlying rotational 
symmetry of the LO state strongly restricts the form of 
the Goldstone-mode Hamiltonian. Under an infinitesi- 
mal rotation of qg, that defines the spontaneously-chosen 
orientation of the pair-density wave, the phase of the LO 
state transforms according to 



9^9 + Sqo-r. 



(68) 



Hence the 9 — —qqu sector of the LO Goldstone-mode 
Hamiltonian must be invariant under this symmetry and 



must therefore be described by a smectic form |l29l4l31 |. 
On the other hand because a rotation of the LO state 
leaves the superconducting phase, (j) unchanged, the (j) 
sector of the Hamiltonian generically does not experience 
any such restriction. We thus expect it to be described 
by a generic anisotropic xy-model form, with the full har- 
monic Goldstone-mode Hamiltonian given by 



niLO 



|(vi.)2 + |(ap)2 



(69) 



Thus, unlike the FF state, the LO state is characterized 
by nonzero, but unequal superfluid stiffnesses, pjr ^ pi. 

Another feature of the LO state is that in addition to 
the primary order parameter, A^o l|64p. it is character- 
ized by a uniform "charge" -4 superconducting and by a 
neutral 2qo-smectic secondary order parameters, 



A(4) 

^SC 



= A 



2|A,j2cos [2q-r + 20], 



(70a) 
(70b) 



where in above we neglected the subdominant contribu- 
tions. As we will see in subsequent sections, these order 
parameters become particularly important when the pri- 
mary order parameter A^o vanishes either due divergent 
fluctuations (as e.g., for T > in two and three dimen- 
sions) or via a disordering transition driven by unbinding 
of topological defects. 



C. Goldstone-mode Hamiltonian 

To support above symmetry-based arguments, we will 
now use the Ginzburg-Landau theory (|50p to explicitly 
derive the Goldstone-mode Hamiltonians for the FF and 
the LO states. 



1. The Fulde-Ferrell Hamiltonian 

To this end we use AF_F(r) inside the T-Lgl ([50|) . but in 
contrast to the earlier mean-field calculation that deter- 
mined the value of A 
Goldstone mode 



qg , focus on the spatially dependent 
r). However, we will neglect the spa- 



tial dependence of the amplitude Aqg , valid in the ordered 
phase, where its deviations from the average condensate 
value are gapped, controlled by a finite susceptibility. 

Working out the gradients of Ai?j;- (r) under these con- 
ditions and using the expression for jpF inside Hgl, we 
find 
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n 



FF 



J|Ao 



(V 



2 J,^2 



(2q.V0+(V(/))2)'' 



(j|A,j2(q2 - ql) + 2^|A,on (2q • V0 + (V0)^ 



ii^(vW + ipii(9||0 + 



(71a) 
(71b) 



In above, we dropped constant pieces, used the FF am- 
plitude Aqj, (computed in the previous section), defined 
a longitudinal derivative d\\ = qo • V, and in going to the 
final form (j71bp chose q according to 



ql 



go, near hc2, 



(72a) 
(72b) 



in order to eliminate the terms linear in the fluctuation- 
current nonlinear form 



5j\\=d\\^+]^q^\V<l>f 



(73) 



as a standard minimization condition for the FF ground 
state current jp f ■ This condition is closely analogous to 
the choice of an order parameter magnitude to eliminate 
terms linear in fluctuations. As expected from the ear- 
lier mean-field analysis, near hc2 (where A^^ is small) this 
corresponds to the choice of g « go- However, as usual 
with a magnitude of an order parameter (here the spon- 
taneous current, Jff), the nonuniversal magnitude of q 
(proportional to Jff) will be modified by fluctuations. 
Ultimately it is determined by the requirement that the 
coefficient of the term linear in Sj^\ vanishes. 

Within this derivation, the Goldstone-mode moduli in 
Hff are given by 



K = 2J|A,J^ 



(74a) 
(74b) 



but more generally are two independent parameters char- 
acterizing the energetics of the single Goldstone mode of 
the FF state. The Hamiltonian form, Hff (valid be- 
yond its weak-coupling microscopic derivation) is fa milia r 
from studies of conventional smectic liquid crvstals |l29l - 
Il3l| . with the rotational invariance encoded in two ways. 
Firstly, to the quadratic order in V(/) it reduces to the 
harmonic form (I6ip . invariant under an infinites- 

imal rotation of the FF current state. Namely, by the 
strict vanishing of the (Vj^i/))^ superfluid stiffness, p-^ = 
(resulting in the "softer" transverse Laplacian energet- 
ics, (V^^c/))^), it exhibits a vanishing energy cost for 



transverse (to qo) current fluctuations, with the stiff- 
ness for the change in current magnitude (along qo) 

controlled by pi. Secondly, Hff is an expansion in 
a fully rotationally-invariant longitudinal current fluc- 
tuation, (5j|| (I73p . whose nonlinearities ensure that it is 
fully rotationa lly invariant even for large reorientations 
in an'l 29l - [T3l| that defines the FF ground state. 

To see the latter we note that under a rotation of qg — 
qoz by an angle a in the z ~ x plane 



qoz 



q = go (-2 cos a + x sin a) 



(75) 



generates a nontrivial, spatially-dependent phase 

0°(r) = z(cosa — 1) + xsina, (76) 

even though the system is clearly in its ground state. 
Simple algebra demonstrates that the fully nonlinear 
form of the longitudinal current dj\\ ensures that it and 
the corresponding energy Hff vanish for 0*'(r), as re- 
quired by the rotational invariance. 

One might question the necessity of keeping nonlinear- 
ities rnHpF (|71bp . As we will see shortly, because of the 
vanishing transverse superfluid stiffness in the FF state, 
the fluctuations in the purely harmonic description H^p 
(1611) are infrared-divergent in three and lower dimensions. 
Consequently these nonlinearities are in fact absolutely 
essential for a well-defined description of such a state. 



2. The Larkin-Ovchinnikov Hamiltonian 

The derivation of the Larkin-Ovchinnikov Goldstone- 
mode Hamiltonian follows a similar route to that for the 
FF state of the previous subsection, with many common 
features, but also some essential qualitative differences 
in the results. To this end, we insert the LO order pa- 
rameter Aio(r) inside Hgl (l50|) . use the earlier mean- 
field values of the amplitude A„q , Eq. ([57)) (that vanishes 
at hc2 and grows as \hc2 — h\° below hc2), ignoring its 
subdominant spatial dependence, and track the resulting 
energetics of two spatially dependent Goldstone modes, 
(j){v) and d{r) = — gou(r). We thereby obtain 



17 



H 



LO 



J|A,J2 f(v2</,+)2 + (2qo . V0+ + {V<t>+ff + (V20_)2 + (2qo • V(^_ - (V0_)2)'' 



2m? 



E 



subdom 
LO ' 



(77a) 
(77b) 

(77c) 



where we dropped the constant and fast oscillating parts 
that average away upon spatial integration of the energy 
density, introduced two phonon fields 



u± = T0±/go, 



(78) 



and chose q = qq in order to eliminate the term linear in 
the nonlinear, rotationally-invariant strain tensor (analog 
of Sj« in dZSl)) 



99 



q- Vu± - i(V'u±) 



(79) 



whose nonlinearities in u± ensure that it is fully rotation- 
ally invariant even for large rotations. We also defined 
the bend K and the compressional B smectic elastic mod- 
uli 



K = 4Jg2|A, 



B = 16Jql\Ag 



(80b) 
(80c) 
(80d) 



and identified the longitudinal, pi 
superfluid stiffnesses given by 



and transverse, 



Ps = 



epqo 



ep 



A, 



90 



^BCS 



epqo 



(81a) 
(81b) 
(81c) 
(81d) 
(81e) 



A nonzero transverse superfluid coupling, (minimized 
by a vanishing supercurrent j oc -I- V(/>_) removes 

the two independent rotational symmetries, orientation- 
ally locking the two incommensurate {u±) smectics. As 
argued above based on symmetry, ([M)) this leads to the 



superconducting phase combination, (f> = ^{(f>+ + </>-) to 
be of a conventional xy- (as opposed to "soft" smectic) 
gradient type, (|77cp . 

Within the superconductor context, this coupling of 
V0± (xy- rather than smectic-like stiffness of the (j> Gold- 
stone mode) is straightforward to understand. The LO 
state effectively corresponds to two FF states, each car- 
rying a supercurrent along q -I- V(/)+ and — q -I- V0_. 
A superfluid stiffness measures the nonzero kinetic en- 
ergy cost of a nonzero total current j. The latter results 
from either a transverse orientational misaligned of V(/)± , 
measured by the p^, or a mismatch between the V0+ 
and V0_ magnitudes (even if directed along ±q, respec- 
tively), measured by the pi. These two distinct LO state 
distortions are schematically illustrated in FiglS) 

As advertised in the introduction we thus find from 
([5T|) that the LO state is a highly anisotropic superfiuid 
(though less so than the FF state, where =0), with 



Ps_ 

pl 



*9o 



^BCS 



4 \ h 



« 1, 



(82) 



c2- 



a ratio that vanishes for h ^ h 

Because (in contrast to the FF state) the p^ is indeed 
nonzero for h < /ic2, the nonlinearities in Vcj) and the 
"soft" term {V^cj))^ (contained in the nf^^'^°'" of (jTTcl) ) 
are subdominant on lo ng scales, and can (and will) there- 
fore be neglected [nil- 

We stress that while the detailed expressions for the 
moduli above are specific to the weak-coupling BCS limit 
near hc2 the general form of Hlo, (|77cI) , including the 
structure of the symmetry-enforced nonlinearities in the 
u (smectic) sector is valid beyond our microscopic deriva- 
tion, and holds throughout the LO phase. Specifically, as 
argued in the previous subsection based on symmetry, in 
contrast to the superconducting, (f) sector, the vanishing 
of the (V_lu)'^ coupling is a refiection of the underly- 
ing rotational invariance of the LO striped state. Thus 
the "soft" transverse elastic form of the smectic sector 
and the form of the nonlinear strain tensor, Uqq (analo- 
gous to (5j|| ([75]) ) are strictly protected by this symmetry. 
As we will see in Sec|V] (and is well-k nown for conven- 
tional smectic liquid crystals |l 29l4l3 1 1 ) . because of this 
harmonic elastic "softness" , in the presence of thermal 
fluctuations these smectic nonlinearities are qualitatively 
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duU^ -duU 




Ji - VjU+-VjU_ 



FIG. 5; An illustration of longitudinal and transverse (to q) 
LO nonzero supercurrent configuration controlled by the cor- 
responding superfluid stiffnesses, pi, p^. Because there is no 
change in the effective period associated with the transverse 
current excitation, it is qualitatively energetically "cheaper", 
scaling compared to the longitudinal excitation that 

scales like A^^. As indicated in EQ. (|82p . the ratio therefore 
vanishes as h approaches the transition into a nonsuperfluid 
phase at hc2- 



important at long length scales for the description of the 
FF and LO states. 



D. Quantum dynamics 

The bosonic Hamiltonian density Hgl ([50]) determines 
the classical thermodynamics and equal-time finite tem- 
perature correlation functions of the FFLO svstem |l54l |. 
that in the ordered state reduces to Hff/lO: in (!71bp 
and ((77c)) . 

The dynamics is controlled by the Lagrangian den- 
sity. Focusing on quantum thermodynamics, the par- 



tition function Z = Trace 



-I3H[A^,A] 



is formulated in 



a standard way as a path-integral over coherent states 
labeled by the c-fields A(t, r) and A*(t, r) 

Z = J[dA*dA]e-^^^''^\ (83) 

with T = it the imaginary time < r < ^9 = l/fc^T 
(fi = 1). The coherent state action S'[A*, A] is given by 



(84) 



with the Lagrangian density given by 

C = A*drA + n[A*,A]. 
Using the density-phase representation 



(85) 



A = (no + (5n+)i/2e*i""+*"^+ + (no + Sn_y^^e"'i-''+"t'- , 

(86) 

with 6n± the fluctuating Cooper-pair density about its 
ground state value hq — |Ag„p, and at low energies 
neglecting the subdominant terms (e.g., fast oscillating 
pieces or spatial dependence of gapped amplitude fields) 
the Lagrangian density reduces to 



C = iSn^dr(t>+ + iSn^drif'- + 'H[5n±, (f) 



±J- 



(87) 



Integrating these Berry's phase terms |l55l Il56| over the 
"massive" (nonzero compressibility, xo) '5"-± fluctuations, 
we obtain 



(88) 



with H given by the H lq in the LO ground state, and FF 
state treated similarly using Hff- For the LO state, this 
analysis then predicts the existence of two anisotropic 
low-frequency modes with dispersions 



w„(k) 



Kk\ 



Bkl, 



(89a) 
(89b) 



that can be read off from the analytical continuation of 
the Lagrangian into real time, t = —it. These modes 
respectively correspond to the zeroth sound (the Bogoli- 
ubov mode as in a conventional superfluid) and smectic 
phonon, unique to the LO state. In cold atomic gases, 
these should in princi ple be me asurable via the Bragg 
spectroscopy techniaue [l57l4l59| 

With the Goldstone-mode Lagrangian in hand, we can 
now calculate the effects of quantum and thermal fluc- 
tuations as w ell as equilibrium correlation and response 
functions [1591. 



IV. LARKIN-OVCHINNIKOV STATE NEAR h^i 

As emphasized earlier the general form of the action 
for the description of the LO state, Eq. ([55)) is expected 
to hold universally throughout the ordered state. How- 
ever, its above derivation and therefore the expressions 
for the associated couplings Egs. dST)) . (|80)) . (|8T)) are lim- 
ited to the weak coupling BCS regime and near the high 
chemical potential imbalance (Zeeman field) normal-to- 
FFLO transition at hc2- 

An estimate of these couplings outside Id throughout 
the FFLO phases can only be done numerically. However, 
in the complementary, low chemical potential imbalance 
(Zeeman field) regime, just above the transition from 
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the fully paired superfluid (BCS-BEC) state to the LO 
state at /ici, qualitative estimates are possible based on 
an analysis of a "dilute gas" of fluctuating ±A domain- 
walls. We carry out this analysis below by focusing on 
the LO state, treating it as a periodic array of fluctuating 
domain- walls in A (r), akin to the lyotropic phases in soft 
condensed matter [l29l Il30l |. 

A. Macro- vs micro-phase separation: stability of 
the LO state 

Implicit in our analysis below is the assumption that 
as the domain-wall surface energy becomes negative [95l- 
losj for h > hci, their interaction remains repulsive, and 
so the domain- walls proliferate continuously as a periodic 
array inside the LO state. Under this assumptions (that 
warrants further study) the domain- wall density ridw and 
the associated species imbalance P oc Udw (~ Qo{h)) 
is then set by a balance between the negative surface 
energy and the domain-wall repulsion, growing continu- 
ously as a function of h — hd according to the Pokrovsky- 
Talapov's commens urate - incommensurate (CI) transi- 
tion phenomenologyjlO^ . This behavior is clearly ex- 
hibited in ld[95i, llOSj through an exact solution and 



where e^^, is the domain- wall surface energy at ft, = 0, 
V{z) is the domain-wall interaction energy per unit of 
area, and rifi^^z) = ~ ^i) {^{^) ^-nd Zi are the 

profile of width £^dw and the position (along z) of the ith 
domain-wall. In above, we used an approximate rela- 
tion between the domain- wall Id density ndw{h) and the 
magnetization density m{h) = n-f — 

m{h) « mNih)^dwndwih), (91) 

where mjqih) is the magnetization density (per d- 
dimensional volume) of the normal state, that is ap- 
proximately nucleated in the zeros of the (locally BCS- 
like) gap function, i.e., on the domain-wall of width £^dw] 
mjs[S.dw is the fermion number imbalance per unit area 
{L'^~^) per domain-wall. In the Pauli (weak imbalance) 
limit, we expect the former to be approximately given 
by kp~^h/hci, through a string of relations: niN^dw ~ 
XpHo ~ inhci/eF)ifiVF/ABcs)ih/hci) ~ kp'^h/hd w 
k'^^ (using n kf, ). This is indeed the case in the 
noninteracting Fermi gas limit, where misi{h) is given as 



bosonization methods, a.nd has been argued to persist 
in higher dimensions |95l - [98| . The CI route for a tran- 
sition to the LO state contrasts sharply with the Lan- 
dau theory 0, 13, [tJ of two independent order parame- 
ters Aq, Ag, that always predicts a first-order BCS-LO 
transition. The latter corresponds to the case of an at- 
tractive domain-wall interaction, that therefore prolifer- 
ate discontinuously above hci, leading to the ubiquitous 
phase separation found in mean- field theory 0, H|. It 
is currently unclear what dimensionless microscopic pa- 
rameter, analogous to Abrikosov's k (disti nguishin g be- 
tween type I and type II superconductors) [93, llOCll | , con- 
trols these two alternatives of the macro-phase separation 
(a first-order transition) and the micro-phase separated 
LO state (a continuous transition out of the gapped SF 
state) [Toll . 



B. SF-LO transition at hd 

The phenomenology of the domain-wall proliferation 
above hd and the associated (fully gapped, singlet) SF to 
LO transition can be captured by a domain-wall energy 
functional, extended to a general dimension d: 



(90a) 
(90b) 



the solution of 

-m-Nih) = -[(A + M^e(A + ft)- (/i-/i)^©(A-/i)], 

(92) 

with the normalized chemical potential jl = /i/e_F de- 
termined in terms of the normalized chemical potential 
difference (Zeeman energy) h = h/ ep by the number den- 
sity equation 

1 = ^[{p + h)iQ{fi + h:) + {fi-h)iQ{p^h)\. (93) 

Above set reduces to the Pauli expression toat ~ Xp^ 
for weak h <ti fj, (xp = 3n/(2ei?)) and « n for large 
h ^ fi. In either case the h dependence of niN^h) is 
weak around hd and therefore can be neglected, along 
with other weak h dependences, such as the soliton width 

Ignoring fluctuations [l60|, the domain- wall interaction 
V{z) sa Vge"!^!/^''"' is expected to be short-ranged with 
the scale set by the soliton width S^dw Thus, for low soli- 
ton density ndwS.dw ^ 1 its strength is a strong function 



E[ndw]/L''- ^ = {Sd^ - hmNih)(,dw)ndw + 7^ / V{z - z')ndw{z)ndw{z'), 

J z ^ J z,z' 

= I mN{h)S,dw{hd~h)ndw + \ I V{z - z')ndw{z)ndw{z'), 

J z ^ J z,z' 
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oi z — z', whose typical value is given by the domain- 
wall spacing itself l/nj^^^. Thus, just above hd a good 
approximation for the above energy density is given by 



£{ndiv) 



h)ndn 
h)ndn 



-^LL LU J ^ — -7 Z 



(94) 



A minimization then gives the domain-wall density, 
ndw{h) and (through Eq. ((9T|) ) the imbalance density 
m{h) as a function of chemical potential difference h: 



ndw{h) 



{h - hci)/ho, for Udwidw > 1, 



f-l J ln[/io/(h-/i,i)] 
^dw ^ 



(95) 



where £,dw 

Vo / {xphci(,o) 

,d-i 



^0 and ho 
Abcs (since 



kp' Abc'5 is the energy lost per fermion per unit of area 
[with separation \/kp\; for a Id exact solution, one in- 
deed has hci — 2Abcs/t^, consistent with this estimate). 

Note that for large imbalance (high domain-wall den- 
sity) the interaction part reduces to a simple quadratic 
dependence on Udw and thus gives the expected linear 
growth of the density with h — hd (second line above). 
In contrast, for low density the interaction is exponen- 
tially weak and domain- walls enter the state as a quickly 
growing (with a divergent slope) function oi h — hd, 
though not as discontinuously as in a first-order tran- 
sition. These limits are illustrated in Figl2l 



1. LO elastic moduli at T = 

For h > hd a Id lattice of domain- walls, i.e., the LO 
state forms. Because of the underlying rotational invari- 
ance its elasticity is generically given by that of a smec- 
tic, Eq. ()77cp . The energy of the deviation Sndw from 
the minimum domain- wall density Udwih), Eq. ()95|) . is a 
quadratic function of Sudw, given by 

e{Sndni) ~ eQ + ^e"{ndw{h)){Sndru)'^, (96a) 
« e„ + ^B{h){d,u)^ (96b) 

where the smectic bulk modulus is given by 

Bih) « nd^ihfe"{nd^ih)), (97a) 
« {h-hd)ndn,ih)k^F''. (97b) 



In above, to deduce B(h) we used a relation between 
the density Udw and the displacement field Sudw = 
—ndwih)dzU, a relation m^r ~ k'^^/$,dw, and went to 
a continuum limit via / ndw{h)dz . . . = X)i • ■ •■ 
expected, B{h) grows strongly just above hd, but 



quickly asymptotes to B{h » hd) ~ Ascskp ^/^q ~ 
uAbcs/^f, hi scale consistent with its form found below 
hc2, Eq. dlQl). 

A detailed estimate of the bend modulus K is more 
difficult in this regime. However, we can take advantage 
of the relation in Eg. (150)) to deduce K via K w B/q^ 
mean-field relation. Thus we find 



K{h) « B{h)/ndUh)\ 



Ascskp ^Co ~ <^Fkp ^, 



for h:$> hd- (98b) 



From these, we deduce the key smectic penetration 
length 



X{h) = ^K{h)/B{h), 
~ l/ndw{h), 
sa ^0, for h > hd- 



2. LO elastic moduli at T > 



(99a) 
(99b) 
(99c) 



Because the microscopic repulsive domain- wall interac- 
tion is short-ranged (exponentially weak at long scales), 
sufficiently close to hd (where solitons are dilute) ther- 
mal fluctuations always qualitatively modify the above 
T = pred ictions. T o understand this we first derive 
the Helfrich |l02l Il6l| interaction between two {d — 1)- 
dimensional fluctuating, curvature dominated (tension- 
less) domain-walls, separated by an average distance z. 
The corresponding energy functional of the instantaneous 
local domain- wall separation w(x) is given by: 



(100) 



where the curvature modulus is related to that of the 
smectic via k = K/udw The average transverse domain- 
wall spacing is determined by the root-mean squared fluc- 
tuations, given by: 



,2\ _ 



T 

K 

T 



Jd-l 



,5-d 



(101a) 
(101b) 



This gives the thermal collision length xt as function of 
separation 



XT 



(- 



l/(5-d) 



2/{5-d) 



(102) 



To deduce the entropic (Helfrich) interaction, we note 
that domain-walls separated by distance z reduce each 
others fluctuation entropy by an amount sq (of order 1) 
per collision. For a pair of domain-walls of linear extent 
Lx there are {L^/ xtY^^ collisions, and thus the entropic 
part of the free energy is raised (relative to the infinite 
separation) by 



5Ft 



Tsc^ ( — 

XT 



d-1 



(103) 
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leading to the Helfrich curvature-controlled interaction 
T" 1 



Vh{z) 



7-4/3 I 

' k1/3 ^2/3' 



with 



/3 
7 



(104a) 

for 2d LO smectic state, 

(104b) 



(105a) 
(105b) 
(105c) 
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Clearly, asymptotically close to the hd transition, 
where the domain-wall array is sufficiently dilute, the 
above thermal fluctuations-induced steric interaction 
Vh{z) always dominates over the short-range microscopic 
interaction. The crossover density nj^ is set by a separa- 



tion 1/nJ^ = 

by 



'■dw 



at which these are comparable, given 



(106) 



The corresponding Zeeman field range is hd < h < hx ~ 
hci + T'^ / lkp~^ ^J^k/^) , with the energy density in this 
regime given by 



mN{h)^dw{hci - h)ndw + —g-nj^^ 



h)ndw + VH{ndw)ndw, 

(107a) 

(107b) 



A minimization then gives ndw{h) 



ndw{h,T) 



for hci < h < hx, 

, (5-d)/2 1/2 



(108a) 



T2/(rf-i) 

for hd <h< hr. (108b) 



As expected it shows that a fluctuation-enhanced Hel- 
frich domain- wall repulsion leads to a significantly slower 
(than the T = 0, Ea. dMjl ) power-law, I/7 increase in the 
domain-wall density and therefore of the species imbal- 
ance P{h). Correspondingly, for this low range of Zeeman 
field h < hx, this enhances the smectic bulk modulus, 
B{h,T) through Ea. ((97b)l . 

For h > hx the microscopic exponential interac- 
tion takes over and the growth of ndwih) and B{h,T) 
crossover to that of the T = result, Eq. ([95)) . 



GOLDSTONE MODES FLUCTUATIONS 



A. Gaussian fluctuations 



As an estimate of the role of Goldstone modes fluctua- 
tions, we first study them at the Gaussian level, namely 
approximate the Goldstone-mode action S[u, 0] at the 
quadratic level, 5'°[w, 0] — 5'"„Ju] + Sg^[(j)] by dropping 
the nonlinearities in Hlo, (I77cp . Combining with ([55)) 
and focusing on the LO state (leaving the straightforward 
extension for the FF state for later) , the LO harmonic ac- 
tion is given by 



oO 



LO 



dT I dzd'^'^r^ 



B 



K 



1 



r 



, (109) 



where (for later mathematical convenience) we general- 
ized the model to d dimensions, with a single LO modu- 
lation ordering axis z = fy along qo and d — 1 space rj^ 

transverse to qo, and introduced n = 2goXo, X — 2xo- 
The harmonic (imaginary time-ordered) correlations of 
the decoupled modes ii(T, r), 0(r, r) can be easily com- 



puted exactly, for 2-point correlation functions giving 

(UOa) 



G„(T,r) = (u(r,r)?/(0,0))o, 
'•^^ d'^k 



R ^ 



(2t:Y kujI + Bkl + Kk\ ' 
G4r,r) = (0(T,r)0(O,O))o, (110b) 
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A Jd 



±1,2 



The averaging above was done with the Euclidean prob- 
ability distribution, e^^^o / Zq, using the harmonic ac- 
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tion above, and w„ = 27rn//3 is the standard Matsub- 
ara frequency. Above, A is the UV cutoff for the Gold- 
stone mode action, set by the inverse coherence length 
1/f ~ Qo- While the full expression above can be com- 
puted asymptotically in terms of special functions or nu- 
merically, it is more revealing to analyze these in spe- 
cial limits of interest. The simplest measure of fluctu- 
ations is given by the root-mean-squared (rms) fluctu- 
ations of these Goldstone modes, given by Gti(0,0) = 
(w2),G4O,O) = (02). 



1. T = quantum fluctuations 



2. T > thermal fluctuations 

At nonzero T the Goldstone modes u and cf) ex- 
hibit classical thermal fluctuations. By the identifica- 
tion of the corresponding sectors with the well-studied 
smectic and anisotropic xy models, we can conveniently 
take advanta ge o f the large body of literature on these 
systems [12^ Il30l Il34j |. However, for completeness we 
will work out some of the key findings. 

Returning to the full expressions for the harmonic cor- 
relation functions, (jllOb[) . (|110cp at nonzero T, we sep- 
arate out the dominant classical Wri=o = contribution 



At zero temperature, the Matsubara summations in 
above expressions reduce to frequency integrals over oj, 
in d dimensions giving 



dojd'^k 



1 




^ dujd'^k 



for d > 1, 



(Ilia) 



(111b) 



A^-i 



Because of the integrand's anisotropics, the details of 
above expressions are sensitive to microscopic ultravio- 
let (UV) and infrared (IR) cutoffs. However, the key 
unambiguous finding above is that in any physical dimen- 
sion of interest here (d > 1) quantum Goldstone mode 
fiuctuations in the LO state are finite, set by the short 
(UV) length scale and therefore are qualitatively unim- 
portant within the ordered LO state. That is, the LO 
state has a nonzero range of stability to quantum fluctu- 
ations, and therefo re at T = is well- approximated by 
its mean-field form [l62j |. Quantum fluctuations will, of 
course, give quantitative corrections to the properties of 
the LO state and will be important near quantum tran- 
sitions to other putative phases. 

Above flndings are consistent with known results for 
mathematically related systems. The action (I109p for 
the d-dimensional T — LO quantum phonon, is iso- 
morphic to a classical Hamiltonian of a generali zed d-l- 1- 
dimensional columnar liquid crvstals |l29l ll3Cll |. charac- 
terized by two spatial "stiff" directions {z and r) and 
d — 1 transverse "soft" rj_ axes. The latter is known (by 
simple power-counting) to exhibit long-range order down 
tod-|-l = 5/2 dimensions, which is consistent with our 
flnding for (w^), above. Similarly, (|110cp is also consistent 
with the well-known property of the xy-model exhibiting 
long-range order down to ci -I- 1 = 2 dimensions. Thus, 
at T = we conclude, as advertised in the Introduction, 
that a d -|- 1 dimensional superfluid smectic, i.e. the LO 
ground state is stable to quantum fluctuations for d > 1. 



d'^k 



T 



L-^ i2nr Bkl -f Kki ' 

^ lngo-^_L, d = 3, 
T 



d'^k 



{"^^Y phi + pikl 



S= — ln9o^_L, d = 2, 



(112a) 
(112b) 
(112c) 

(112d) 



where we neglected the subdominant quantum contri- 
bution (worked out above), defined a constant Ca = 
Sd/{2nY = 27r''/V[(27r)''r(d/2)], with Sa a surface area 
of a d-dimensional sphere, and introduced an infrared 
cutoff by considering a system of finite extent L± x L^, 
with Lz the length of the system along the ordering (z) 
axis and L± transverse to z. Unless it has a huge aspect 
ratio, such that Lz ~ » L±, any large system 

{L_i,Lz >> A) wiU have XL^ < L^.- 

The key observation here is that the smectic phonons 
exhibit fluctuations that diverge, growing logarithmically 
in 3d and linearly in 2d with system size L±. The su- 
perconductor phase fluc tuations 6 exhi bit well-known 2d 
logarithmic divergences [ill, Il39l Il4lj . Because (aside 
from the superfluid stiffness anisotropy) these are the 
same as in an ordinary superconductor and in the bulk 
3d case will be finite, indicating a long-range off-diagonal 
order, we will focus on the smectic phonon fluctuations 
special to the LO (and FF) states. 

The expression for the root-mean squared phonon fluc- 
tuations in (|112b|) leads the emergence of important 
crossover length scales , ^2 , related by 



6 



(e./?^^B)^/^ 



(113a) 
(113b) 



that characterize the flnite-temperature LO state. These 
are defined as scales L±, at which phonon fluctuations 
are large, comparable to LO period a = 27r/qo- Namely, 
setting 



,2\T 



(114) 
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in Eq. (ITT2bl) we find 



BK 



_K_ 

iira'^^/BK/T , 



(115a) 



where in the second form of the above expressions 
we took the simplest approximation for the smectic 
anisotropy length A = \/ KjB to be A = a ^ l/ff n, and 
introduced an order 1 Lindemann constant c[163[, that 
depends on the somewhat arbitrary definition of "large" 
phonon rms fluctuations. 

The thermal connected correlation function of LO 
phonons 



where erf(a;) is the Error function. 

Above finding of the divergence of smectic phonon fluc- 
tuations at nonzero temperature have immediate drastic 
implications for the properties of the LO (and FF) states. 
As emphasized in the Introduction, the most important 
of these is that the thermal average of the Landau's LO 
order parameters ((64)) vanishes in thermodynamic limit 



Cu(y^,z) - ([u(ri,z)-u(O,0)]')o 



(116) 



is also straightforward ly w orked out, in 3d giving the 
logarithmic Caille form jl34| 
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(2^)3 Kq\_ + Bql ' 
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In — - -Ei , 

a J 2 V4A|z| 
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where Ei(a;) is the exponential-integral function. As in- 
dicated in the last form, in the asymptotic limits of 
r±_ >• \/Az and rj_ <C VXz above 3d correlation func- 
tion reduces to logarithmic growth with r± and z, re- 
spectively. 

In 2d we instead have 1421 




(Aio(r))o = 2A,„(e^'^cos[qo-r + (?])o, 



2A, 



qa'^ 



COS (qo • r), 



2A5o(i_L)cos (qo • r). 



(119) 



with the thermally suppressed order parameter ampli- 
tude given by 



(117a) 
(117b) 



0, for L± -> oo, 



(120b) 



where (j)rms = (0^)o- In above we used results for the 
phonon and phase fluctuations, (jll2b|) . (jll2d|) and de- 
fined the Caille exponent 



V = 



(121) 



We also neglected the subdominant quantum phonon 
fluctuations and included all finite (quantum and ther- 
mal) superconducting phase fluctuations inside the 
nonzero Debye- Waller factor, e~5<^'rms. Thus, in qualita- 
tive contrast to its mean-field description, at long scales 
(longer than f_L,z) a LO state is characterized by a uni- 
form superconducting order parameter and density. 

To further characterize a finite-temperature LO state, 
we compute the Cooper-pair momentum distribution 
function, 
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(AIo(r)ALo(r'))e*-""-^'^-""', 



C*(r)g-igJC„(r)g»(k-9„5).r 



(122a) 
(122b) 

(122c) 



r 



where ^ is the system's volume and we generalized the 
LO order parameter to include high harmonics g„. To 
lowest order we neglected the subdominant coupling be- 
tween the superfluid phase and the LO phonons, ignored 
the nonlinear effects, approximating u and (j) as Gaussian 
fields. The asymptotic form of n^'^ can then be readily 
obtained analytically. As we saw above, in 3d the super- 
fluid phase, (j) exhibits only finite fluctuations about an 
ordered state, with 



Cf{r » 



T 



(123) 



where we took the coherence length ^ as the short-scale 
cutoff. These phase fluctuations then simply give a finite 
Debye- Waller factor suppression of the momentum dis- 
tribution function amplitude, corresponding to the usual 
fluctuation-driven condensate depletion, with a fraction 
of atoms pushed out of the qo condensate. 

In contrast, the 3d LO phonon fluctuations diverge 
logarithmically, (|117b|) . strongly modifying rik from its 
mean-field ^-function form. We thus find to exhibit a 
power-law peak around the ordering wavevector qo (and 
its harmonics, q„), reminiscent of (l-l-l)d Lutti nger liq- 
uids and two-dimensional crystals 136l - ll38L Il4ll | 



nqo 



\2 — n^ri ' 



for d - 



(124) 



where for simplicity we specialized to k = fc^z. The 
form-factor amplitude is approximately given by k, 



An additional characterization of the LO state is 
through a structure function, S'(q), a Fourier transform 
of the density correlation function, that in 3d is given by 



'S'^°(q) = y dV((5p(r)<5p(0))e-''i'-, (125a) 

« 1 dV(|ALo(r)nALo(0)|')e-^'^•^ (125b) 
« ^|A,J4 jr(e'29'.("o-".))oe-''(i-2g„£).r^ ^^25c) 
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where we approximated phase and phonon fluctuations 
by Gaussian st atist ics (in 3d valid up to weak logarith- 
mic corrections jl 3 ij ) and replaced atomic density fluctu- 
ations by (twice) the LO condensate density. The lat- 
ter neglects a contribution from the imbalanced atoms, 
without qualitatively modifying the result (since atomic 
density is locked to the condensate density) and is fur- 
thermore quantitatively subdominant for h <C hc2^ where 
the imbalance is low. As for nt, we find that the loga- 
rithmically divergent 3d phonon fluctuations lead to a 
structure function, with highly anisotropic (gz ~ ^l/-^) 
quasi-Bragg peaks (see Fig i2T|) replacing the true Bragg 
peaks characteristic of the mean-field long-range periodic 
order. These p redictions are a reflectio n, of the well- 
known jl29lll3d l and experimentally tested |l35j | behavior 
of conventional smectic liquid crystals. 



In two dimensions, the LO order is even more strongly 
suppressed thermal fluctuations. The linear growth of 
the 2d phonon fluctuations leads to exponentially short- 
ranged correlations in the LO order parameter and in the 
density. Because it is the "soft" smectic Goldstone mode 
that is responsible for these interesting properties they 



• qper t 



are necessarily also shared by the FF state 



As discussed in the Introduction, another fascinating 
consequence of the thermal vanishing of (A^o), (jl20bp is 
that the leading nonzero Landau order parameter char- 
acterizing the LO state is the translationally-invariant 
"charge" -4 (4- atom pairing) superconducting order pa- 
rameter, Asc, introduced in (I70ap . Thus in the pres- 
ence of thermal fluctuations the LO phase corresponds 
to an exotic state in which the off-diagonal order is ex- 
hibited by pairs of Cooper pairs, i.e., a bound quartet of 
atom s, ra ther than by the conventional 2-atom Cooper 
pairs |ll9l |. In 2d and 3d this higher order pairing is 
driven by arbitrarily low-T fluctuation, rather than by 
a fine-tuned attractive interaction between Cooper pairs, 
and therefore has no mean-field description. We will dis- 
cuss these and other fluctuation-induced phases as well 
as transitions between them in SeclVIIl 
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B. Nonlinear elasticity: beyond Gaussian 
fluctuations 

1. Perturbation theory 

As is clear from the derivation of the previous subsec- 
tion, the restoration of the translational symmetry (a uni- 
form LO state with a vanishing (Aio(r)), etc.) by ther- 
mal fluctuations is a robust prediction of the quadratic 
theory, that cannot be overturned by the left-out nonlin- 
earities. However, the asymptotic form of the correlation 
functions computed within the harmonic approximation 
only extends out to the nonlinear length scales £,±.z, be- 
yond which the divergently large LO phonon fluctuations 
invalidate the neglect of the nonlinear phonon operators 

Hnonlinoar = -^B {d,u)iV uf + {V u)\ (126) 

These will necessarily qualitatively modify predictions 
(fTT7b)) . (fTT8B . (ITM)) . and (fHS)) on scales longer than the 
crossover scales that we compute next. 



6B above) become comparable to the bare microscopic 
values. The divergence of this correction as L± oo sig- 
nals the breakdown of the conventional harmonic elastic 
theory on length scales longer than a crossover scale 




, d — 3 



which we define here as the value of L± at which 
\dB{^^^)\ = B. Within the approximation of the smectic 
screening length X — a, these nonlinear crossover lengths 
reduce to the phonon disordering lengths (|115p . (jll3bp . 
defined by a Lindemann-like criterion. Clearly, on scales 
longer than the perturbative contributions of nonlin- 
earities diverge and therefore cannot be neglected. Their 
contribution are thus expected to qualitatively modify 
the harmonic predictions of the previous subsection. 



2. Renormalization group analysis in d = 3 — e dimensions 




FIG. 6: Feynman graph that renormalizes the elastic moduli 
K, B of the LO superfluid. 

To see this, we use a perturbative expansion in the non- 
linear operators (|126|) to assess the size of their contribu- 
tion to e.g., the free energy. Following a standard field- 
theoretic analysis these can be accounted for as correc- 
tions to the compressional B and bend K elastic moduli, 
with the leading contribution to SB, summarized graph- 
ically in Figini and given by 
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(127b) 



In above, we neglected the subdominant contribution 
from quantum fluctuations, i.e., used thermal equal-time 
correlator, G'„(q), focused on d < 3 (which allowed 
us to drop the uv-cutoff (A) dependent part that van- 
ishes for A oo), and cutoff the divergent contribution 
of the long wavelength modes via the infra-red cutoff 
q_L > 1/Lj_ by considering a system of a finite extent 
L±. Clearly the anharmonicity become important when 
the fluctuation corrections to the elastic constants (e.g.. 



To describe the physics beyond the crossover scales, 
- i.e., to make sense of the infra-red divergent per- 
turbation theory found in Eq |127b1 - requires a renormal- 
ization group analysis. This was first performed in the 
context of conventional liquid crystals and Lifshitz points 
in a seminal work by Grinstein and Pelcovits (GP)[l3l|. 
For completeness, we complement CP's treatment with 
Wilson's momentum-shell renormalization group (RG) 
analysis, extending it to an arbitrary dimension d, so 
as to co nnec t to the behavior in 2d, that has an exact 
solution 132 [ . 

To this end we integrate (perturbatively in Hnoniinoar) 
short-scale Goldstone modes in an infinitesimal cylindri- 
cal shell of wavevectors, Ae^^^ < q± < A and — oo < 
^2 < oo {S£ ^ 1 is infinitesimal). The leading perturba- 
tive momentum-shell coarse-graining contributions come 
from terms found in direct perturbation theory above, 
but with the system size divergences controlled by the 
infinitesimal momentum shell. The thermodynamic av- 
erages can then be equivalently carried out with an effec- 
tive coarse-grained Hamiltonian of the same form (|77cp . 
but with all the couplings infinitesimally corrected by the 
momentum shell. For smectic moduli B and K this gives 
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where dimensionless coupling is given by 



g = C,_iAl-r ( 
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(130a) 
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and in the second form we approximated g by its value 
in 3d. Because it is only the smectic nonhnearities that 
are qualitatively important, all other stiffnesses do not 
experience corrections that accumulate at long scales. 
Eas. (|129l) show that B is softened and K is stiffened by 
the nonhnearities in the presence of thermal fluctuations, 
making the system effectively more isotropic. 

For convenience we then rescale the lengths and the re- 
maining long wavelength part of the fields u"^ (r) accord- 
ing to r_L = r'j^e^^ z = z'e'^^'^ and M<(r) = e"^'"^u'(r'), 
so as to restore the ultraviolet cutoff Aj_e~^^ back up 
to A_L . The underlying rotational invariance insures that 
the graphical corrections preserve the rotationally invari- 
ant strain operator {pzU — ^(V_lm)^), renormalizing it as 
a whole. It is therefore convenient (but not necessary) to 
choose the dimensional rescaling that also preserves this 
form. It is easy to see that this choice leads to 







(131) 



The leading (one-loop) changes to the effective coarse- 
grained and rescaled action can then be summarized by 
differential RG flows 

'^^^^^ ^ {d+-i--iu^\g{i))B{l) , (132a) 
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dl 
dK{i) 



1 



id-l-uj+-g{i))K{i) . (132b) 

From these we readily obtain the flow of the dimension- 
less coupling g{£) 



(133) 



whose flow for d < 3 away from the g = Gaussian fixed 
point encodes the long-scale divergences found in the di- 
rect perturbation theory above. As summarized in FigjT] 
for d < 3 the flow terminates at a nonzero fixed-point 
coupling g^, = (with e = 3 — d), that determines the 
nontrivial long-scale behavior of the s yste m (see below). 
As with treatments of critical points [l6J|, but here ex- 
tending over the whole LO phase, the RG procedure is 
quantitatively justified by the proximity to d = 3, i.e., 
smallness of e. 



lation functions on scales beyond Ci'^. Namely, apply- 
ing above coarse-graining RG analysis to a computation 
of correlation functions allows us to relate a correlation 
function at long length scales of interest to us (that, be- 
cause of infrared divergences is impossible to compute 
via a direct perturbation theory) to that at short scales, 
evaluated with coarse-grained couplings, B{£), K{£),. . . . 
In contrast to the former, the latter is readily computed 
via a perturbation theory, that, because of shortness of 
the length scale is convergent. The result of this match- 
ing calculation to lowest order gives correlation functions 
from an effective Gaussian theory 



G„(r = 0,k) 



T 



BCk)k^,+K(k)kl' 



(134) 



with moduli i3(k) and i^(k) that are singularly 
wavevector-dependent, latter determined by the solutions 
B{e) a nd K{e) of the RG flow equations and 
(|132bp with initial conditions the microscopic values B 
and K, e.g., as given by the BCS predictions, Eq. (jSO]) . 

2d analysis: In d = 2, at long scales g{i) flows to a 
nontrivial infrared stable fixed point g^ = 32/5, and the 
matching analysis predicts correlation functions charac- 
terized by anisotropic wavevector-dependent moduli 

K{k) = if (fcieD"""/K(fc.^r/(^±eD'^), (135a) 



B{\.) = B{k,_i^^y^ fsik^er/ik^er)'), (135b) 

Thus, on scales longer than these qualitatively 

modify the real-space correlation function asymptotics 
of the harmonic analysis in the previous subsection. In 
Eqs. p35p the universal anomalous exponents are given 
by 



1 4 

8 
4 



for d = 2 , 
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FIG. 7: Renormalization group flow for a LO state in d < 3- 
dimensions, illustrating that at low T it is a "critical phase" 
displaying universal power-law phenomenology, controlled by 
a nontrivial infrared stable fixed point. 

We can now use a standard matching calculation to 
determine the long-scale asymptotic form of the corre- 



for d = 2 , 



(136b) 



determining the z — rj_ anisotropy exponent via (jl34p to 
be 



C = 2-(?7s+?7k)/2, 
7 

" 5' 



(137a) 
(137b) 



as expected reduced by thermal fluctuations down from 
its harmonic value of 2. The k_L — kz dependence of 
i3(k),_ft'(k) is determined by universal scaling functions, 
fsix), fK{x) that we will not compute here. The under- 
lying rotational invariance (special to a LO state realized 
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in an isotropic trap) gives an e xact relation between the 
two anomalous r/g x exponents [l65j 



This translates into an equal-time LO order parameter 
correlations given by 



3-d 
1 



2 

VB_ 

2 



r]K , for d = 2, 



(138a) 
(138b) 



which is obviously satisfied by the anomalous exponents, 
Eqs. (ll36bp . ()136a|) . computed here to first order in e = 
3-(iil65i. 

Thus, as advertised in the Introduction, we find that 
a finite temperature 2d LO state is highly nontrivial and 
qualitatively distinct from its mean-field perfectly peri- 
odic form. In addition to a vanishing LO order parameter 
and associated fluctuation-restored translational symme- 
try, it is characterized by a universal nonlocal length- 
scale dependent moduli, Eq. (I135p . Consequently its 
Goldstone mode theory and the associated correlations 
are not describable by a local field theory, that is an ana- 
lytic expansion in local field operators. Instead, in 2d, on 
length scales beyond thermal fluctuations and cor- 
relations of the LO state are controlled by a nontrivial 
fixed point, characterized by universal anomalous expo- 
nents riK^B and scaling functions fB,K{x) defined above. 

Above we obtained this nontrivial structure from an 
RG analysis and estimated these exponents within a con- 
trolled but approximate e-expansion. Remarkably, in 2d 
an exact solution of th is problem was discovered by Gol- 
ubovic and Wang 132| . It predicts an anomalous phe- 



nomenology in a qualitatively agreement with the RG 
predictions above, and gives exact exponents 



J2d 

Vb 
Vk 

A2d 



1/2, 
1/2, 
3/2. 



(139a) 
(139b) 
(139c) 



n(z,r^-0) - (Alo(r)Aio(O)), 



Loy-Z^^oy^)/, (141a) 
3-"i('"^)'''cos((jo2), (141b) 



(ci a nonuniversal constant) as discovered in the co ntex t 
of conventional smectics by Grinstein and Pelcovits (l3l| . 
Although these 3d anomalous effects are less dramatic 
and likely to be difficult to observe in practice, theoreti- 
cally they are quite significant as they represent a quali- 
tative breakdown of the mean-field and harmonic descrip- 
tions, that respectively ignore interactions and thermal 
fiuctuations. 

We conclude this section by noting that all of the above 
analysis is predicated the validity of the purely elastic 
model, Eq. (|77cl) . that neglects topological defects, such 
as vortices and dislocations. If these unbind (a s the y 
undoubtedly do in 2d at any nonzero temperature [m^), 
then our above prediction only hold on scales shorter 
than the separation between these defects. 



VI. TOPOLOGICAL DEFECTS IN A 
LARKIN-OVCHINNIKOV STATE 



3d analysis: In = 3, the nonlinear coupling g{i) is 
marginally irrelevant, flowing to at long scales. De- 
spite this, the marginal flow to the Gaussian fixed point 
is sufficiently slow (logar ithm ic in lengths) that (as usual 
at a marginal dimension [l64l | ) its power-law in £ depen- 
dence leads to a universal, a syrn ptotically exact logarith- 



mic wavevector dependence [131 



if(ki,fc, = 0) ^ X|l + -^ln(l/fcia)p/5 

647r 



B(ki =0,fc,) - B\l 



5^ 
Mn 
5.9 
1287r 



ln(A/fc,a^)| 



2n|-4/5 



(140a) 
(140b) 



We now turn to the discussion of the topological de- 
fects, followed in the subsequent section by an analy- 
sis of phases and transitions accessible by their unbind- 
ing. As discussed in Sec lIIII the LO superfluid is distin- 
guished by two independent order parameter components 
^±qoj correspoi iding to ±qo finite center of mass mo- 
mentum pairing |l52j . These complex order parameters, 
A±qg = |Ag(,|e"^± in turn lead to two independent phase 



Goldstone modes 0± (or equivalently 4',9 ~ ~qou, 
controlled by a long-scale Hamiltonian, that at harmonic 
level is given by 



= / dzd''-^rj_ 



E ( T^9.M' + ^(Vi0.)^ ) + ^(Vj.<^+ + V^0_) 



4%^ 



8 



(142a) 
(142b) 
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with the couphngs given in Eqs. (l5D)) . (|8ip and in obtaining 
Eq. (I142p we dropped the subdominant higher gradient 
{K) term in (p. The analysis of this Hamiltonian in the 
absence of topological defects was discussed in the previ- 
ous section. We now use it to understand the energetics 
of defects beyond that "spin-wave" approximation. 

As in an ordinary superfluid, because (f)± are compact 
phase fields and (j)± + 2tt are physically identified), 
in addition to their smooth configurations, there are vor- 
tex topological excitations, corresponding to nonsingle- 
valued configurations of 4>±{r). These are defined by two 
corresponding integer-valued closed line integral enclos- 
ing a vortex line 

(j) d£-^(t)±=2TTn±, (143) 

that we collectively designate by a two component 
integer-valued vector Ny — {n^,n^), with n± G Z. 
These integer vector defects, iV„ are associated wit h the 
fundamental group Hi of the torus [/(I) (g) t/(ll. jl66l | 
that characterizes the low-energy manifold of Goldstone 
modes of the LO state. In this respect the LO superfluid 
has similarities to other t7(l) Uj l) systems, such as 
easy-plane spinor-1 condensates[167] and two-gap super- 
conductors, e.g., MgB9 il6& |. 

In a differential form, the line defects are equivalently 
encoded as 

V X V(t)± = m± , (144) 
with vortex line topological "charge" density given by 

m±(r) = 27r^ / nl^U{s.,)S^ir - r^{s^))ds^ , (145) 

i 

where Si parameterizes the i'th vortex line (or loop), 
ri{si) gives its positional conformation, ti{si) is the local 
unit tangent, and vortex "charges" rij. are independent 
of Si, since the charge of a given line is constant along 
the defect. Furthermore, 

V • m(r) = (146) 

enforces the condition that vortex lines cannot end in the 
bulk of the sample; they must either form closed loops or 
extend entirely through the system. 

A. Vortices and dislocations 

As with the Goldstone modes in ([M]). where it was 
more convenient to work with the more physical sum and 
difference modes, 0(r) = i(0+ -!-(/)_), 6{r) = i((/)+ — (/>_), 
Eqs. (|65| . (|142bp we consider topological def ects associ- 
ated with singularities in (f>{r) and 0(r),'l66| defined by 



(j) d£- = 27rn^, (147a) 

jdl-Ve = 27md- (147b) 

Given the definitions of </>, 9, the corresponding vortex (v) 
and dislocation (d) defects "charges" , are related to 
the integer-valued n± according to 

riy = i(n+ + n_), (148a) 

Ud = i(n+ - n_), (148b) 

and therefore admit half-integer and integer topological 
"charges" , collectively designated by a vortex flavor vec- 
tor A/" = [ny^n^). The minimal values of n+ ~ ±1, 
?i„ = ±1, lead to four fundamental defects (8 counting 
the overall ± signs of each) 

a/; (±1,0) ^ 7V„,i, = (±1,±1), (149a) 

A/d = (0,±l) o Ar,^_, = (±l,Tl), (149b) 

A/'^^d = (±l/2,±l/2) o A?„.o = (±1,0), (149c) 

A/^;_d = (±l/2,Tl/2) ^ A^o,. - (0,±1), (149d) 

where the flrst two defects, A/'„,A/'d are ordinary 27r- 
vortices in 0+ and in , that respectively correspond 
to a 27r- vortex in superimposed on defect-free LO lay- 
ers (FiglSl) and an integer a-dislocation in the periodic 
LO smectic state. 

As illustrated in Fig.[5]for a 2d state, the latter of these 
two corresponds to a LO state with an edge dislocation 
of a smectic layer, or equivalently, two missing adjacent 
± domain-walls, with no additional phase winding in (j). 
In 2d the LO order parameter in the presence of these 
point defects is given by 

A^O = 2e*'^(^'^)cos(go^), forA; = (l,0), (150a) 
Aio = 2cos[qaz + ed{y,z)], for A/^ - (0, 1), (150b) 

where Lp{y,z) — tan^^(z/j/) is the azimuthal coordi- 
nate angle creating the singular vortex, and Od{y,z) = 
ip{y,x) + 0o{y,x) the corresponding dislocation angle 
with a nonsingular part 9o{y,x) accounting for smectic 
anisotropic elasticity. As can be seen from their form in 
terms of the = (n+,n_) description, these two seem- 
ingly simpler integer defects in Eas. (|149aj) . (jl49bp are ac- 
tually composites of the fundamental J\fv,Q — (1,0) and 
7Vo,t) = (0, 1) defects. 

The novel A/^„_d ^ (±1/2, ±1/2), At^-d = 
(±1/2,=f1/2) defects are half-integer vortex-dislocation 
composites, where a ±7r- vortex is bound to a half of single 
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FIG. 8: A 27r-vortex defect in a LO state. The 2d arrows indi- 
cate the corresponding complex condensate LO wavefunction, 
A5|,o(r), with orientations characterizing the phase <^(r). 
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FIG. 9: An integer a-dislocation defect in the LO layered 
structure, with no defects in its superfluid phase. 



FIG. 10: Two complementary forms illustrating a half-integer 
vortex-dislocation (7r,a/2) defect in the LO state. The 2d 
arrows indicate the corresponding complex condensate LO 
wavefunction, A^Q''(r), with orientations characterizing the 
phase 0(r). 



in the presence of such a half-integer defect is given by 

A^o** = 2e^^(^'-)/2cos[qoz + ^(y,z)/2], (151) 
for a/; = (1/2,1/2). 

In contrast to a conventional uniform superfluid or a 
standard smectic density wave, here the product form 
of the pair-density wave LO order parameter (I64p allows 
this composite half-integer defect. Namely, although the 
superfluid and the smectic density wave order-parameter 
components each change by a minus sign (wind by a 
phase of tt), the product LO order- parameter A^o, 
Eq. ()64p remains single-valued in the presence of such 
(±1/2, ±1/2), (±1/2, =Fl/2) fractional composite defects. 



B. Energetics 

domain- wall, an a/2-dislocation, illustrated in FigHO] In 

terms of the two coupled phase fields, 0^., these half- The energy cost of above described vortices and dislo- 

integer defects correspond to an integer vovtex in one, cations can be comp uted with straightforward extensions 
but not both 0± phases. In 2d the LO order parameter of standard analvses |l30l Il42| . 
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1. 2tv -vortex: AJ^ = (1,0) 

The 27r-vortex in the LO superfluid is obtained via a 
solution of the Euler-Lagrange equation for 0(r) 



(152) 



under the condition of a 27r-vortex line singularity 

VxV0" = 2n I i{s)S^{r-r{s))ds , (153) 



located at r(s) with a unit tangent t(s). The smectic 
layers remained undistorted at long scales, u = (equiv- 
alently = (j)-)- The solution displays a qualitatively 
standard form, but with distortions associated with a 
large anisotropy ratio p\\/p±, that, as we demonstrated 
dramatically diverges at the upper-critical field hc2- This 
is reflected in the anisotropy of the mass current flow 
around a vortex line running within the LO (smectic, 
xy-) layers, (a mass flow with velocity component along 
the smectic layer normal, qo = z), illustrated in FigfTT] 
For concreteness we look at two straight vortex line con- 
figurations. 

2tt -vortex line ||qo- 

For a straight vortex line running along qo (z-axis) 
above equations reduce to 

Vi0" = , (154a) 
V X = 2TTzS^{r±_), (154b) 

which gives ^''(r) — (p{x,y), where ip{x,y) = taii~^{y/x) 
is the azimuthal coordinate angle within a smectic layer 
(xy-plane). The corresponding superfluid "velocity" 
= V^" for such a 27r- vortex line directed along qo = ^ 
(flowing in the xy _L-plane) is isotropic, given by a stan- 
dard ip/r± form 

^l^^hy,x,0). (155) 



Integrating the kinetic energy density, Eq. (ll42bD m a 
system of dimensions L± x L± x Lz we readily find the 
energy to be given by the familiar 3d (linear-logarithmic) 
form, 



El ^TrpjLzlniL^/a), 



(156) 



diverging linearly with vortex line's length, and loga- 
rithmically with systems in-plane extent L±. 
2'K-vortex line _L qo: 

In contrast, a vortex line directed along the smectic 
layers (taken here to be along the x-axis, with flow con- 
fined to the anisotropic y — z plane) characterized by 

V X V^" = 2T:x5{y)5{z) (157) 



is described 



tan 




by (/)(r) = ip{ 
and a velocity field 



PsV; 



PsPi 



p^z^ 



(0,-z,y), 



(158) 



illustrated in FigfTT] 

As expected from the superfluid-stiffness anisotropy 
near hc2, where pjr /pi <C 1, the superflow around a vor- 
tex falls off much more slowly across the LO layers (along 
z) than within them (along rj^). This leads to (near hc2 a 
highly) anisotropic elliptic-like superflow with the major 
axis oriented along the LO layers normal (along qo) and 
the velocity flow that is not everywhere tangent to the 
equipotential hows, as illustrated in FigfTT] 




FIG. 11: Anisotropic superfluid flow around a 27r- vortex near 
hc2, predicted to be characterized by anisotropic superfluid 
ratio pl/p't that diverges at hc2- 



Substituting this flow inside the Hamiltonian (I142bl) . 
the corresponding vortex energy is then readily shown to 
be given by 



E^=nJplpiL^\n {pUl + p^Ll^^^. 



(159) 



In 2d it reduces to a familiar logarithmic form by replac- 
ing the Lj_ factor by the sample thickness. 



2. a-dislocatton: Nd = (0, 1) 

An integer dislocation defect in the LO superfluid is 
similarly obtained via a solution of the Euler-Lagrange 
equation for u = —O/q^, obtained by minimizing 
Ea ll42b[ with a condition of cj) — Q (equivalently 0+ = 



Bdlud - K^\ud = , 
and the singularity and continuity conditions 



V X vd 



0, 



(160) 

(161a) 
(161b) 
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with the singular strain 

Vrf = Vud , (162) 
and the dislocation "charge" density 

md(r)=a^ J nl^ii{si)S^{r - r,{si))dsi . (163) 

i 

As for the vortex, Si parameterizes the z'th dislocation 
loop, ri(si) is the conformation of that loop, ti(si) is 
its local unit tangent, and is the number of excess 
layers associated with the dislocation. Note that is 
independent of Si, since the charge of a given dislocation 
line is constant along the line defect. 

The dislocation solution is most straightforwardly ob- 
tained in Fourier space, where equations reduce to: 



and the general solution is given by 



Vd = 5 h iqx(q) 



(164a) 
(164b) 



(165) 



with x(q) a smooth elastic distortion around the dis- 
location line obtained via the Euler-Lagrange equation, 
(|160p . and is given by 



g.(l-AM) , 
^ q y 



(166) 



The superfluid phase remains uniform at long scales, = 
(equivalently = —</)_). In above we defined the 
inverse of the LO phonon propagator 



-psm _ 2 I \2 4 
i q = + ^ 9j. 



(167) 



with A = y^K/B, a length characterizing the LO domain- 
wall deformations. Inserting above expression for 
into the original elastic Hamiltonian Ea ll42bl gives the 
Coulomb gas-like Hamiltonian that determines the en- 
ergy of an arbitrary dislocation configuration 



Hd = 



1 



Kql 



q 



/^|(q)+£e4 



77if(q)m^(-q), (168) 



where Pj^(q) = (5^ ^^I'tlj' ll\ iii-(xy) plane trans- 

verse projection operator, with _L indicating projection 
onto the smectic layers. In above we also added dislo- 
cation core energy determined by its short-scale (a — 
27r/(7o) configuration, at which above continuum analy- 
sis no longer applies. A simplest estimate (up to possibly 
nontrivial dimensionless factors set by A/a, that can only 
be determined through a microscopic calculation) of this 
core energy density £c is given by K. 

In 3d, two limiting dislocation line configurations are 
the screw and edge dislocations. 



edge 



'■(r) 
(r) 



= az(5^(rj_), 
— axS{y)S{z), 



(169a) 
(169b) 



running perpendicular to and along smectic layers, re- 
spectively. Within above harmonic continuum approxi- 
mation, the long-scale elastic contribution (first term in 
(|168p ') to the energy of a screw dislocation vanishes iden- 
tically (due to P^{q)), physically because of smectic's 
soft in-plane curvature elasticity and due to the absence 
of long-scale compressional distortion in a screw disloca- 
tion, as can be readily verified. Thus its energy is de- 
termined by short-scale core energetics inaccessibility to 
harmonic elastic analysis. 

For a single edge dislocation along x (|169b|) m'*(q) = 
ax2Tr6{qx), giving for the elastic part of the energy 



i3d—edge 

d.el " 



/ ^(l-'Z^/gi)2vr%.), (170a) 
q q 



Ka 



dq-ydqz Kq: 



V 



(27r)2 q^^+XX'' 



(170b) 
(170c) 



where the factor of Lj_ came from the regularized iden- 
tity 27r(5((73;)(5((72,) = L±d{qx)- Although the first "elastic" 
contribution in Eq. (jl68|) no longer vanishes, the above 
energy is dominated by the short cutoff scale, a and there- 
fore in principle indistinguishable from the core energy. 
Thus, based on this analysis we can only provide an es- 
timate of the 3d dislocation energies. 



E 



3d— screw 



E' 



<3d—edge 



KL 
Ka 



(171a) 

L_L ~ {BKf/'^aL^ w KL^, (171b) 



diverging linearly with its length. 

Corr espo ndingly, as first demonstrated by Toner and 
Nelson [l42|, in two dimensions the energy of a smectic 
edge dislocation is constant, i.e., system length indepen- 
dent, given by 

Ef w ay/BK^K. (172) 

As we will see in the next section, this observation has 
crucial implicatio ns for the instability of a finite temper- 
ature 2d LO statefiii. 



3. half-integer vortex-dislocations defect: Nv-d = (1/2, 1/2) 

As discussed aboveflT^, on general grounds a LO 
(but not a FF) superconductor allows novel half-integer 
defects that are a composite of a tt- vortex and a/2- 
dislocation, illustrated in FigHOl The form and the en- 
ergy associated with this defect requires a solution of 
Euler-Lagrange equations for (l)± 



<i 



-p^Vi^^ = 0, 
(173) 
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with singularity conditions 



where 



V X V0+ = 27:xd{y)d{z), (174a) 

V X V<?!)_ = 0, (174b) 



for concreteness taken to be a straight vortex line in (/>_)_ 
running parallel to LO layers, along the x axis, and (f>- — 
0. 

Solving above equations gives the defect's shape from 
which the energy is readily computed 



E. 



7r-a/2 ~ 4^-^ 



.^In [{plLl+piLl)'/ya]+K 



Because the energy scales as a square of the defect's topo- 
logical charge n^^d, as expected, the above energy of the 
half vortex-dislocation defect is one-quarter of the sum 
of the energies of the unit vortex and dislocation. 



C. Dual description: coupled xy-smectic models 

In characterizing phases in terms of topological defects, 
it is often useful to have a continuum field-theoretic de- 
scription of defects and of the corresponding transitions 
associated with their unbinding. Such a description, that 
complements the LO Goldstone model, (I142p is obtained 
through a duality transformations that we develop here. 
In 3d the analysis is somewhat complicated and results 
in two coupled U{1) gauge theories, whos e der ivation and 
analysis we leave to a future publication [l5l|. Here, for 
simplicity we will exclusively focus on two dimensions. 

To this end, starting with Hlqi (|142bp . allowing for 
singular configurations of (j) (vortices) and u (disloca- 
tions) discussed in the previous subsection, and integrat- 
ing out the smooth parts of Goldstone modes, we obtain a 
Coulomb gas description of a finite density of topological 
defects 



jjxy— sm 
^CG 



^ q 



^C{,V I 
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^q,d| 



(176) 



r^y — 



psm 

q 



pUps + qzV Ps/p. 



(177a) 
(177b) 



In above, mq^t,,mq_d are the Fourier transforms of the 
vortex and dislocation densities m„(r) = 2Tin'^.5^{Y — 
r^),md(r) =^,a<,52(r-rf), with n^f ^ \{n+±n^X 
and n^. g Z the independent elementary integer defects, 
defined in the previous subsection. In above we also in- 
(175) troduced the vortex and dislocation core energy £'"^^ to 



account for short scale core physics not accounted for 
by the above continuum description. Although by time- 
reversal symmetry the core energies of the elementary de- 
fects are identical, those of the vortex and dislocation 
composites will generically be distinct with 7^ Ef. 

To obtain a continuum description of these integer- 
valued fields, it is convenient to decouple interaction of 
these defects using a Hubbard-Stratonovich transforma- 
tion by introducing a pair of dual real fields (correspond- 
ing potentials) 9 with the Hamiltonian 



i7[0,6l,n„,<j] = Ha[4>,d]+i / ( 0(r)m^(r) + qo6'(r)TOd(r) 

(178) 



where 



yxy 

piPs 



T^sm J2 

t3_Jl\a |2 



(179) 



The thermodynamics is characterized by a partition 
function 



z = I [d4>de] n E e-"^^''^<^ 

v,d v,d 



(180a) 
(180b) 



r 



where we chose to measure all energies in units of ksT. can be readily done in the dilute defect limit, correspond- 
The summation over vortex and dislocation charges n'^f ing to large core energies Ec- In this limit the summa- 
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tions can be limited to nine lowest order terms, corre- 
sponding to no vortices present, Nq_q, one vortex of ei- 
ther species and sign at a site, N±i,o, -^o,±i, or two vor- 
tices of distinct species sitting at a single site, N±i ±i. 
Equivalently, in terms of M vortex-dislocation nomencla- 
ture, these respectively correspond to 



This gives 



A^: (0,0), ±(7r,a/2), ±(7r. 



-a/2), ±(27r,0), ±(0,a). 

(181) 



[d4>de] 



.-Ho 



1 + 2^2 (e-3^"-3^' cos [7r{4> + 9)] 



+e~^^ cos(27r(/)) -f- e cos(27r. 



e~i^c~iEi cos [t:{4>-9)] 



(182) 



r 



All other vortex configurations (for example two defects 
of either species sitting on distinct sites) are not included 
as they correspond to higher order contributions, sup- 



pressed at low fugacity. After re-exponentiating, above 
leads to a generalized sine-Gordon model for the coupled 
(h, 9 fields, with the dual Hamiltonian 



H = - 
2 



q 



psm^2 



ff7r,a/2 cos(7r0) cos(7r6i) -I- 527i-,o cos(27r^) 



+ go,a cos{2tt9) (183) 



r 



where the couplings are (7,r,a/2 



527r,0 



and go,a 



4e- 



and we converted the sum 



over ri into an integral over r via a lattice constant a. 

In the opposite limit of a small core energy, E'^!''^ the 
defect density is high and the summation over the inte- 
ger charges nf can be carried out utilizing the Poisson 
summation formula, that leads to a replacement of the 
e.g., cos(27r0) potentials above by the Villain potential, 
defined by: 



-l3Vv 



^e'^^^4>-E.n^ ^ (184a) 

n 



The latter can be approximated by a single harmonic 
with an effective coupling gl^ = l/{2E^-'^a?). 

The resulting generalized sine-Gordon model is con- 
venient for analyzing the effects of defects on the LO 
state, particularly for the computation of their screening 
on long scales, unbinding, and for the analysis of the re- 
sulting disordered state. From the form p83p it is clear 
that (aside from an inconsequential anisotropy) the dual 
vortex sector described by has a standard sine-Gordon 



form. In contrast, the dual dislocation sector, described 
by 9 is qualitatively modified by the highly nonlocal and 
qualitatively anisotropic smectic kernel. 



We next turn to a detailed discussion of implication of 
these defects for possible phases and transitions associ- 
ated with the LO state. 



VII. PHASES AND PHASE TRANSITIONS: 
LARKIN-OVCHINNIKOV LIQUID CRYSTALS 



We now turn to a discussion and characterization of 
a class of phases emerging from the Larkin-Ovchinnikov 
state by a partial disordering of it. Because these phases 
partially break spatial symmetries, we will refer to them 
collectively as LO liquid crystals. We first discuss gen- 
eral symmetry based possibilities and then explore their 
concrete realization in terms of fluctuation-driven prolif- 
eration of topological defects. 
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Landau's order parameters and spontaneously 
broken symmetries 



A conventional characterization of phases is through 
the Landau's classification, where phases are distin- 
guished by non-vanishing order parameters and corre- 
sponding symmetries that they break. From this prospec- 
tive a zero-temperature LO (smectic pair-density wave) 
superfluid {SFsm) breaks three symmetries: the transla- 
tions Tqg along qo (z), the rotations i?qp about the axis 
transverse to qo (of the full 3d Euclidean group E{3)), 
and the U{1) symmetry associated with atom conserva- 
tion. The corresponding non-vanishing order parameter, 
that transforms non-trivially under these group of sym- 
metries is the LO pair-density wave 

By considering all possible basic combinations of spon- 
taneously broken subset of these symmetries, we uncover 
five additional atomic liquid crysta l pha ses, that are de- 
scendants of the smectic LO state jl70l |. To enumerate 
these systematically we begin with the zero-temperature 
LO state, where all three above symmetries are broken 
and partially restore them by progressively disordering 
the state. Restoring the translational symmetry Tq^, 
while keeping U{1) and Rq^ leads to a state with ori- 
entational and off-diagonal orders, that is a nematic su- 
perfiuid, SF^. By analogy with the conventional (non- 
superfluid) nematics the resulting state can be character- 
ized by a complex traceless symmetric second-rank tensor 

Qij • 

Subsequently disordering the orientational order and 
thereby also restoring the rotational symmetry i?q„ 
leads to an isotropic superfluid, SFj, that exhibits a 
finite species imbalance. Symmetry-wise the resulting 
state is isomorphic to the polarized superfluid, SFm 
predicted d, nil and observed to appear at a nonzero im- 
balance on the BEG side of the BCS-BEC crossover. In 
contrast, (as a descendant of the LO state expected to 
be stabilized by Fermi surfaces imbalance) here the SFj 
state is realized in an imbalanced superfluid on the BCS 
side, something that has been searched for dating back to 
Sarmafzil, but has not been possible within mean-field 
treatments, that instead predict an instability to phase 
separationjlliO,!!!]. 

As summarized by Tables IIIIII the additional three 
phases, smectic, nematic and isotropic Fermi liquids, 
FLsm, FLn, FLj are the non-superfiuid counterparts of 
the three discussed above, obtained by first restoring 
the U{1) symmetry by disordering the off-diagonal long- 
range order. The fully disordered FLj state is simply the 
normal state of the polarized Fermi gas, albeit strongly 
interacting. Together these intermediate fluctuation- 
induced phases (along with a number of other possible 
ones that we discuss below) naturally interpolate between 
the fully gapped singlet (homogeneously and isotropic) 
BCS superconductor at zero imbalance and low temper- 
ature, and the normal polarized Fermi liquid at large 
imbalance and/or high temperature. 



FLsm 

|U(1) 

SFs^ 



FLn 

Tu(i) 

SFn 



FLi 
|U(1) 
SFi 



TABLE I: Five pliases that naturally emerge as disordered 
descendants of the LO (superfluid smectic, SFsm) state. 
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FLi 
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FLn 
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FLsm 
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SFi 


X 


V 


V 


SFm 


X 
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X 


SFsn, 


X 


X 


X 



TABLE IL A summary of LO liquid crystal Fermi-liquid (FL) 
and superfluid (SF) phases, and corresponding order parame- 
ters and broken symmetries, indicated by X's. Unbroken sym- 
metries are marked by check marks. The subscripts /, N, Sm 
respectively indicate the Isotropic, Nematic and Smectic or- 
ders. 



Larkin-Ovchinnikov liquid crystals via 
topological defects unbinding 



A complementary way to characterize phases and 
phase transitions between them is in terms of topological 
defects proliferation[l4lj. For conventional phases and 
transitions (e.g., 3d xy or Ising models) this is simply a 
complementary description that is sometimes convenient. 
However, for topological phases, where Landau order pa- 
rameter is unavailable or just insufficient to distinguish 
two phases, this defect-proliferation approach is indis- 
pensable and therefore superior to the order-parameter 
Landau's "soft-spin" description. One prominent and fa- 
miliar example is the description of the low-temperature 
quasi- long-range ordered phase of the 2d xy model and its 
disordering. In this case, both the low- ("ordered") and 
high-temperature (disordered) phases exhibit a vanish- 
ing xy order parameter, and are respectively only topo- 
logically (not symmetry) distinguished by bound and 
unbound vortex defects and by the associated behavior 
of the correl ation functions (power-law and exponential, 
respectively) |l4l| . 

With this in mind, we characterized the LO and its 
descendant states in terms of a proliferation of the four 
types of topological defects N = (0, a), (27r, 0), (tt, ±a/2), 
introduced and analyzed in Sec I VII 

As discussed above, at T = the LO phase is char- 
acterized by long-range off-diagonal and translational or- 
der and is thereby distinguished by a nonzero LO order 
parameter A^o. In this ground state all of the above 
topological defects are absent, confined into topologically 
neutral pairs. 
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1. 3d phases and transitions 

Although, (as we exphcitly demonstrated in Sec|Vl at 
nonzero temperature the LO order parameter vanishes in 
3d (and in 2d) , the 3d LO phase is distinguished from its 
more disordered descendants by the absence of unbound 
topological defects, in direct analogy with the quasi-long- 
range ordered state of the 2d xy model. Thus, in this 
same sense at nonzero T the 3d LO smectic is a topologi- 
cally ordered (but phonon, elastically disordered) phase. 

Upon increasing temperature, or decreasing the stiff- 
nesses (e.g., by tuning the strength of the Feshbach- 
resonant interactions, or by adjusting the fermionic 
species imbalance), one or more of the four topologi- 



where for simplicity we ignored anisotropics in p\ . Based 
on this energetics one may be tempted to conclude that 
in this limit (unless preempted by a first-order transition) 
it is the integer dislocation loop defects that proliferate 
first and the LO smectic preferentially disorders into a 
nematic superfluid, SF m- However, in contrast to the 
2d KT mechanism [l4l| . the 3d disordering transitions 
take place when the relevant stiffness, renormalized by 
quantum and thermal fluctuations is continuously driven 
to zero at the transition, or takes place at a finite (rather 
than a vanishing) fugacity. For a thermal transition this 
roughly corresponds to a transition temperature set by 

the corresponding stiffnesses, — \J pip^ and K, B. 

Lacking a detailed quantitative theory of such 3d tran- 
sitions we can only construct a qualitative phase dia- 
gram, that we display in Fig|31 It summarizes all basic 
phases that naturally appear upon disordering of the LO 
smectic superfluid by unbinding four fundamental type 
of topological defects discussed above. One important 
physical input is the observation that increasing the Zee- 
man energy, h (or equivalently, the species imbalance, 
m = — ni) toward the normal Fermi liquid state, FLj 
at hc2, predominantly leads to a suppression of the su- 
perfluid stiffness and therefore to a destruction of the SF 
order. Conversely, a reduction of the Zeeman field (and 
species imbalance) toward a conventional isotropic and 
homogeneous superfiuid SFj at hd, primarily leads to a 
reduction of the elastic moduli of the smectic pair-density 
wave by increasing its period 1/qo and thereby weaken- 
ing the interaction between the LO domain-walls (see 
Eq. (j97bp ). Thus lowering of h is expected to predomi- 
nantly suppress, i.e., melt the positional smectic order. 

Thus, starting with the LO SFsm state and decreasing 
h leads to the unbinding of the integer dislocations (0, a), 
and a transition to an orientationally-ordered, i.e., a ne- 



cal defects will proliferate, thereby leading to a tran- 
sition of the LO smectic SFsm (PDW) into one of 
its descendants. The actual sequence of defects un- 
binding is determined by the relative energetics, given 
by Eas. (|159p . (11711) . (11751) . Lacking a reliable calculation 
(throughout the P—l/kpa phase diagram) of depen- 
dences of the stiffnesses on the experimentally tunable 
parameters (T, imbalance P, interactions 1/kpa, and the 
number of atoms), the phase diagram can only be quali- 
tatively mapped out in terms of the three effective stiff- 
nesses fiT, pi'"*". In the thermodynamic limit (L_l,z — >■ oo, 
not necessarily relevant to atomic traps; see below), the 
relative defects energetics is quite unambiguous: 



(185) 

I 

matic "charge" -4 superfluid, SF^. The later "charge" -4 
feature of SF^ naturally appears as the remaining sec- 
ondary order parameter Ai'^ — A|^q once the LO posi- 
tional order A lq is destroyed by unbinding of integer dis- 
location loops. General arguments predict a subleading 
singularity ^ {h — h[^f_sm)^~°' of A'fJ{h) at the continu- 
ous transition h^q-Sm between SFj^ and SFsm- Since in 
contrast, the "charge" -2 SF nematic order vanishes in the 
LO state, a direct transition to it from the LO state can 
generically only proceed through a first-order transition. 



Conversely, we expect that increasing h starting with 
the LO SFsm will lead to a suppression of , a prolifera- 
tion (27r, 0), and a transition to a 2(7-smectic Fermi liquid, 
FL^g^, a non-superfluid periodic state with a wavevec- 
tor that is twice the LO state. Alternatively, as sug- 
gested by the energetics in Eq. (|185p . it maybe that the 
lower-energy half-vortex dislocation defects (tt, a/2) (or 
the (tt, —a/2), but not both) unbind flrst, in which case 
a transition to a nematic Fermi liquid, FL*j^ (with the 
restored translational and U(l) charge symmetries) takes 
place. The resulting state is qualitatively distinct from 
the more conventional nematic (orientationally ordered) 
FLn phase in which both (tt, a/2) and (tt, —a/2) are pro- 
liferated. Both are also distinct from the nematically 
ordered FL*jy state, in which only integer dislocations, 
(0, a) and integer vortices, {2tt, 0) are unbound. One can 
envision a number of other states and phase transitions 
at low h by further considering the disordering of the ne- 
matic superfluid, SF^ by unbinding various patterns of 
disclinations and 7r-vortices. We leave a more detailed 
analysis of these to a future study. 



I 

Eto,a) - KL « iJ^--^^/^) ^ In i + |l « El,^^,^ ^ p,L In L, for L^,, ^ L ^ oo, 
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FIG. 12: A flowchart of superfluid (SF) and nonsuperfluid 
(FL) phases, exhibiting smectic (Sm) and nematic (TV) con- 
ventional orders as well as topological orders (indicated by * 
and **), induced by a proliferation of various combination of 
topological defects, (0, a), (27r,0), and (7r,±a/2). 

2. 2d phases and transitions 

The nature of the phase diagram changes qualitatively 
in two dime nsion s at nonzero temperature. Based on the 
work of Ref ll42l (in the context of conventional smectic 
liquid crystals) and from the analysis in Sec lVI Bl that 
shows a finite energy cost of a 2d smectic dislocation, we 
conclude that at any nonzero T in 2d, the dislocations 
proliferate, thereby destroying the LO {SFsm) phase. It 
is replaced by a homogeneous, but quasi-orientationally- 
ordered, "charge" -4 superfluid nematic, SF^. Upon ro- 
tation this superfluid will display 1/4 "charge" vortices 



C. Fractionalized phases and topological order 

In above discussion we argued for the existence of at 
least three topologically distinct Fermi liquid phases that 



ifVcj) ■ dl = TTh/2m), that, because of its nematic order 
we expect to form a uniaxially distorted hexagonal lat- 
tice. Upon changing T, h, kpa, the nematic superfluid 
can then undergo further disordering transitions toward a 
polarized Fermi liquid and an isotropic (e.g., BCS) super- 
fluid. In particular, we expect two Kosterlitz-Thouless 
transitions associated with the loss of 2d superfluid and 
orientational (nematic) quasi-long-range orders. Because 
Aq and therefore ^/5^p| ~ vanishes strongly near 

hc2, we expect the superfluid KT transition to precede 
the nematic-isotropic one. 

I 

naturally emerge from disordering of the LO (SFsm) 
phase by unbinding different combinations of allowed de- 
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fects. Because, as demonstrated above the conventional 
vortex (27r, 0) and a conventional dislocation (0,a) are 
composites of the fundamental defects (tt, ±a/2) the non- 
superfluid states FL*j^, FL*j^ and their isotropic cousins 
FLj, FL*j* (in which disclinatio ns ar e also unbound) are 
expected to be "fractionalized" [l7lj . topologically dis- 
tinct from their conventional Fermi liquid analogs, where 
(tt, ±a/2) are also unbound. 

These novel phases are analogous to the putative 
phase-disordered fractionalized states obtained by un- 
binding double (hc/e) vortices, studied extensively 
by Sachdev, and by Balents, Senthil, Fisher, and 
collaborators 17lMl73j | in the context of high tempera- 



ture superconductors. The resulting nonsuperfluid phase 
is distinguished from a conventional Fermi liquid by a 
gapped "vison" , a Z2 defect that is a remnant of the fun- 
damental hc/2e vortex after the composite hc/e (double) 
vortices proliferate. These states are also characterized 
by "fractionalized" charge-neutral spin-1/2 bosonic and 
charge-e spinless fermionic excitations. 

The states FL*j^, FL*^ also bare a close relation to the 
collecti ve mode fractionalization discussed by Sachdev, 
et al. jl74l4l76} in the context of quantum paramag- 



netic phases, emerging from disordering a collinear spin- 
density wave. As with the {U{1) ® U{\)) /Z2 LO state, 
where the order parameter is a product of the superfluid 
and smectic order parameters. Eg. ([Ml) , there too the or- 
der parameter is of (5*2 U{1)) /Z2 product form, encod- 
ing spatial modulation of the spin density, and therefore 
admits half-integer (vison-like) defects. 

As we have seen in the previous subsection, a charac- 
terization of both the conventional and the fractionalized 
(topologically ordered but otherwise disordered) phases 
can be faithfully formulated directly in terms of distinct 
patterns of proliferation of the four types of vortex and 
dislocation defects, listed in Eq. (|149p . It can also equiv- 
alently be done in terms of a dual sine-Gordon (and in 3d 
U{1) gauge-theory) model H , Eq. (I183|) . We now present 
a complementary effective Ising gauge theory description, 
that can sometimes be convenient. In the absence of 
an underlying rotational invariance (that otherwise has 
been our focus throughout) and for simplicity ignoring 
anisotropy of the striped LO state, all the phases and 
transitions can be captured by a Euclidean action (in 
space-time x = (r, r)) 



J 



^x ^x' 



<7x,x' C0s{9s - Ox') - Kcy'^Yl. '^S,x' 

{x,x') □ 



(186) 



r 



where the Ising gauge field, ag^g' couples to the U{1) ® 
U{1) bosonic (rotor) matter fields e"^^,e'®^, and its 
nonzero Z2 flux through a plaquette encodes the presence 
of a half-integer (tt, a/2) defect. The gauge field as.S' en- 
codes the local Ising redundancy of splitting the LO or- 
der parameter (l64|) into a "charge"-2 boson, b]. — e^"*"', 
that creates a zero-momentum Cooper-pair (diatomic 
molecule) and a neutral boson, pj ,. = e^^^' , that cre- 
ates a density wave at a LO wavevector q. Formally, the 
gauge field a can be introduced as a Hubbard-Stratonovic 
field that decouples the LO order parameter Aq — bpq 
into its charge and density wave parts. The variety of 
phases and transitions between them are summarized in 
the phase diagram, FigUand a flow-chart, FigfT^ 

In above formulation, the LO superfluid {SFsm) is a 
state at large t^^tg and arbitrary K„, in which both b 
and pq are Bose-condensed, and cr is gapped through a 
Higg's mechanism. In this state, gapped 27r vortices in 
<f> and 6 respectively correspond to the (27r, 0) supercon- 
ducting vortex and the (0, a) integer dislocation, that we 
discussed in the previous subsection. 

For small and large te 2it vortices in proliferate, 
driving b normal and restoring the U^{1) (atom conserva- 
tion) symmetry, while keeping pq condensed. A large K^; 
forces a vanishing Ising flux with cr = 1 , that corresponds 
to a gapped vison. The resulting nonsuperfluid state 



is deconfined in a sense that that it exhibits a gapped 
bosonic b excitation carrying an Ising charge and thereby 
acquiring a phase tt upon encircling a vison. In our ear- 
lier notation this is the nonsuperfluid periodic state we 
dubbed FL^g^, in which (27r, 0) vortices have prolifer- 
ated, but dislocations remain bound. Lowering drives 
visons gapless, corresponding to a proliferation of the 
(tt, a/2) fractional defects that induces a transition to 
the homogeneous but orientationally ordered (nematic) 
nonsuperfluid state, FLm- We note that condensation of 
(tt, a/2) in the presence of unbound (27r, 0) defects auto- 
matically also leads to a proliferation of (tt, —a/2) defects 
and therefore (aside from the nematic conventional order) 
the resulting FL^ is fully disordered. 

In the opposite regime of large and small tg 2n vor- 
tices in 6 proliferate, driving pq normal and restoring the 
Ug{l) translational symmetry, while keeping b condensed. 
For large K^^ vison remains gapped and the resulting su- 
perfluid homogeneous state exhibits gapped density exci- 
tations Pq, carrying an Ising charge and nontrivial statis- 
tics with the vison. The resulting state is the charge-4 
nematically ordered superfluid, SF^. Upon reducing K„ 
visons proliferate, driving a transition to FLjq through 
this alternate route. 

In contrast to lowering K^r, the transitions out of the 
FL^g^ and SFf^ states at large K„ can be driven by 
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respectively lowering and te and thereby unbinding 
the second set of integer defects, (0,a) and (27r, 0), re- 
spectively. Since visons remain gapped, the resulting 
nonsuperfluid nematic state is the topologically ordered 
FL*j^i, qualitatively distinct from FLj^. The dcconfine- 
ment transition FL*j^-FL n is then driven by lowering 
through a condensation of visons and is in the inverted 
Ising universality class. 

A naive attempt at a generalization of the above Ising 
gauge theory action, (|186p to the rotationally invariant 
smectic form suggests a replacement of the tg operator 
by a gauge-invariant lattice Laplacian. However, on gen- 
eral grounds, without fine-tuning such Ising lattice form 
appears to preclude a fully rotationally-invariant formu- 
lation necessary for a fuUy-rotationally invariant LO (su- 
perfiuid smectic) state. 

Finally, we note that above orientationally ordered 
state can further disorder into isotropic states by prolifer- 
ation of disclinations. We leave a more detailed study of 
these phases and the corresponding transitions to future 
research. 



VIII. FERMIONS 

So far all of our discussion following the defining micro- 
scopic model in SecHIl has been confined to the bosonic 
sector of the Larkin-Ovchinnikov state. The resulting 
low-energy quantum thermodynamics is encoded in the 
Hamiltonian, Hlq, Ea. (j77cp (and the Lagrangian Clo, 
Eg . (1551) ) for the superfluid phase (j) and smectic phonon 
u Goldstone modes. In contrast to the fully gapped su- 
perconductors where this is sufficient at low energies, the 
gapless nature of the LO state also requires the inclusion 



of gapless fermionic excitations for a complete descrip- 
tion of the state. While a detailed analysis of these is 
beyond the scope of the present manuscript, below we 
comment on a few a key features of the fermionic sector 
of the LO state. 



A. Gapless fermionic excitations in the 
Larkin-Ovchinnikov state near hc2 

As we have seen in our discussion of the Goldstone 
modes, there are two complementary descriptions of the 
FFLO states, respectively valid near hc2 and hd- Just 
below hc2 at the continuous FL-FFLO phase the pair- 
ing order parameter, A(r) is small with a weak sinu- 
soidal modulation, and the momentum-space descrip- 
tion is most appropriate. The corresponding mean-field 
many-body wavefunction for e.g., the FF state is given 
by Ea. (|20l) . It is of the BCS form with a range of k over 
which the energy of the corresponding Bogoliubov quasi- 
particles is driven negative by the Zeeman field. These 
therefore form a Fermi sea with all the phenomenology 
associates with the gapless fermionic excitations at the 
Fermi surface. 

For a more general FFLO state defined by a set of 
reciprocal lattice vectors q„, no exact analytical solution 
of the BdG Hamiltonian is available. The periodic A(r) 
couples all particle (atoms) and hole states connected by 
the reciprocal lattice q„, leading to an anomalous A(r)- 
dependent band structure of Bogoliubov quasi-particles. 

However, some progress can be made for a purely sinu- 
soidal LO state with only ±q reciprocal lattice vectors. 
As discussed earlier a BCS-like mean-field Hamiltonian 
for such a state is given by 



k.a k 
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The difference from the FF BCS Hamiltonian, analyzed 
exactly in SeclH] is in the appearance of both q and — q 
Fourier components of A(r), taken to be equal for the 
simplest cosine LO form. These couple each creation op- 
erator c^j^^ to two fermionic atom annihilation operators 

^k+qt '^k-qf ^ith each in turn coupling to another 
set of two operators at different k's, thereby generat- 
ing an infinite-dimensional space that needs to be diago- 
nalized. Focusing on the three-dimensional subspace for 
each k, the Hamiltonian can be rewritten in the following 



Bogoliubov~-de Gennes form: 

/a+qt Aq 

= E*k^BdG*k + E^-k;, (188) 



with the three component generalization of the Nambu 
spinor given by 



-k+qt 

^-ki 



(189) 
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and HBdG the corresponding Bogoliubov-de Gennes 
Hamiltonian matrix. The nontrivial eigenstates encode 
that the excitation in the LO state is a hnear combina- 
tion of a hole (— k, |), an atom (k + q, t)i and an atom 
(k-q,t). 

Although it appears that the problem reduces to a di- 
agonalization in this three-dimensional subspace, in fact 
(in contrast to the single Fourier component FF state, 
where -ffedc is block-diagonal, here) there is a coupling 
between a bottom component of a spinor and the top 
component of the spinor \E'k'=k-2q, 



(190) 




k-3qt 



This leads to the aforementioned infinite-dimensional 
space to diagonalize, corresponding to the band structure 
of the Bogoliubov quasi-particles, analogous to a system 
with a diagonal periodic potential. 

Despite these complications, for a large enough q (such 
that /2m ^ Ag, well satisfied near hc2 where q sa 1/^) 
an approximate treatment is possible because the cou- 
pling is dominated by the degenerate particle-hole states 
near the Fermi surface, which reduces the problem to 
only a pair of states for every value of k. For positive 
(negative) k the pair is the top (bottom) two compo- 
nents of the spinor ^k, leading to a Cooper-pair with q 
(— q) center of mass momentum. 

Diagonalizing the BdG Hamiltonian then leads to two 
sets of the Bogoliubov quasi-particle operators dko-q, 

Q^k(7— q 



Oktq 

d^ 

dkt-q 
-H-q 



Wk^qCj^+a^ + ^'k,qC_k+a^> (191a) 

-<,qCk+5t + "k.qC^k+^V (191b) 
"k -qCi^_q^ + Wk -qcl^ a^, (191c) 

-<-qCk_at + <-qC-k-a^, (191d) 



with the corresponding four branches of excitation spec- 
trum, E'kcrqi 
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(192d) 



— ^ ~ + ^^'^ coherence factors Uk, Vk, 



approximately given by the FF expressions in Eq. (|15p . 
In the first [second] pair of equations, Eqs. (|191"a|) . (|191b[) 
[Eqs. (ll91c| . (ll91dl) ]. the particle- hole hybridization is via 
Aqg^i '' [Age-*i ""], as for the FF state. 



As demonstrated for the FF state in Sec.(jTT]), it is clear 
from the spectra i^ko-iq, Eq. (|192p . that the LO state ex- 
hibits k regions of both gapped and gapless fermionic ex- 
citations, with closing of the gap driven by a combination 
of the Zeeman (imbalance) energy h and the Doppler shift 
k • q. For pairing driven by A^e*'* '" (A^e^^'i '") the gap- 
less states appear at the minimum of the E^^q (i^kt-q) 
located in a wedge —9m<d<0m around k || — q 
(k II q), where for positive h, -Bktq (-E-k-f-q) is driven neg- 
ative. The LO ground state thus takes the form given in 
Eq. ipPl) . exhibiting Fermi pockets of Bogoliubov quasi- 
particles, with a Fermi surface of gapless excitations de- 
fined by £^kyo.±q = 0, as illustrated in FiglTSl 





FIG. 13: An illustration of Fermi pockets (full curve) of 
the gapless Bogoliubov quasi-particles characteristic of the 
Larkin-Ovchinnikov ground state. The periodic array of 
domain- walls in /\lo{z), the associated wavevector q, and 
the Fermi-surface of the underlying normal state (dashed cir- 
cle) are also indicated. 

In the complementary wedge of k around k - q = 0, i.e., 
running along the LO stripes given by the nodes of the 
LO order parameter, the superconducting gap remains 
finite, growing to its maximum value Aq — /i > (sup- 
pressed by h) as 
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The quadratic k± dispersion along the stripes is char- 
acterized by an enhanced effective mass me// ~ 
mep/ABcs > m. 

While above analysis gives some understanding of the 
LO states, it is approximate and does not address the 
detailed nature of a more generic LO-like states. These 
may be understood by self-consistently diagonalizing the 
BdG equations in real space. An estimate of the mean- 
field BCS-LO and LO-FL phase boundaries requires a 
computation of the total energy as was done in Ref|3 
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for the FF state. A complete solution for the LO state 
requires a numerical analysis [95l - [97| . 



For a single "- to +" domain-wall the solution can be 
readily found (l77l Il78l | and exhibits two normalizable 
Andreev zero-energy bound states, 



B. Gap less fermionic excitations in the 
Larkin-Ovchinnikov state near hd 

The above calculation of the LO quasi-particle spec- 
trum near hc2 can be complemented by a treatment 
valid near hd. At low h, the LO order parameter 
A^oiz) is given by a periodic (period 27r/ o) array of 
±Ao ~ ±Abcs domain- walls of width f |95l - l98l |. rather 
than a single harmonic. 

The form of the quasi-particle excitations is deter- 
mined by a two-component Nambu wavefunction ip — 
{ip^^x., z),ipi{:x., z)) satisfying the BdG equation 



(194) 



where p = V and we took the order parameter to be real 
(characteristic of the mean-field LO state) stripe domain- 
walls lying in the x plane, with the normal along z, and 
approximated the chemical potential by a Fermi energy 
ep, valid deep in the BCS regime. In a matrix form, the 
BdG equation is given by 




P 

2m 



A(z)ct, 



V'(x,z) = E^{:>i,z). (195) 



For a striped LO state we can utilize translational invari- 
ance along the stripes, x, taking z) = 'ipii_^{z)e''^^'^, 
with V'ki {z) satisfying 

2 



2m 



A{z)a^ 



V'k_L (^) 



(196) 



with eF{k±) = ep — k\/2m the effective Id Fermi energy. 

The low-energy spectrum is determined by the excita- 
tions near the Fermi energy 



V'k,(z)«0^Jz)e^'^-^ + ^kJz)e 



-ikp-z 



(197) 



with kp/2m = eF{k±) and the envelope wavefunctions 
0^ (z) satisfying 



±vfPz + -z e_F 

Zm 



A(z)(7, 



(-) 



(198) 

For ep{kj_) sufficiently large, so that we can linearize 
around the Fermi points ±fci? by neglecting the quadratic 
correction to the dispersion, Hp2 



■^u,, we obtain 



[{TivFdz - ep) cFz + lA(z)ax\ 0^ (z) = Ek^c/)^^ (z) 



(199) 



^± {z) ^ A{ \ ] fo Mz')dz' ^ ^200a) 




(200b) 



where ^ = vp/Aq is the effective coherence length (set- 
ting the width of the bound state) and in the second 
line we approximated the "- to domain-wall by a step 
function A(z) — AoSgn(2;) valid for z ^ ^. For the op- 
posite sign "- to domain-wall, A(z) = — Aosgn(z) the 
normalizable Andreev zero modes are given by identical 
expressions, but with the reversed sign in front of the i 
in the second component of the Nambu spinor. 

The minimum excitation gap is determined by the non- 
linear correction to the free spectrum, the perturbation 

Hp2 = l^'^zi that we incorporate through a degener- 
ate perturbation theory in the two-component degener- 
ate subspace of Andreev zero modes, (f'^^iz). A sim- 
ple computation gives the off-diagonal matrix elements 



modes to = ±Ek 



2me ' 

, with 



Eh 



^0 



2m^2 Aep{kA_) 



splitting the Andreev zero 



for fc2^/2m < cf, (201a) 



^0 



4(eF-A:f/2TO) 4eF 



A^ 

2ep 



2m 



(201b) 



The maximum of the gap can be estimated through a 
variational solution of the Eq. (jl98p for the nondegenerate 
case of vp — Ip — 0, giving 



rpmax 

gap 



^o(Ao/eF)l/^ 
Ao(Ao/eF)'/' 



(202a) 
(202b) 



for the width of the state localized on the domain-wall 
and the maximum of the gap. 

For a periodic array of domain-walls with period a = 
27r/(7o (the LO state), the coupling between the An- 
dreev states localized on each domain wall splits their 
energy into bands separated by gaps and dispersing with 
—ir/a < kz < Tr/a. The eigenstate for the bottom of 
the band, kz = kp can be obtained by generalizing the 
single domain- wall zero-energy state, Eq. (|200l) to a peri- 
odic array Alo{z)- The solution for the bottom of the 
Andreev band is simply a sum of two periodic arrays of 
states, localized on "- to (at z = 0) and on to -" 
(at z = a/2) domain-wall arrays, respectively 
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4A^^ L 



^-ih lo '^i'0{z')dz' 



Ia/2 ^LO(z')dz' 



(203) 



For a nonzero chemical potential difference h (Zeeman 
energy) to impose the spin imbalance, the spin up ex- 
citation spectrum (I201bl) simply uniformly shifts down 
by h. h exceeding the minimum gap = Af^/Aep, 

Eq. (|20ip . generically induces Fermi pockets (Fig|T3| of 
the Bogoliubov quasi-particles cente red around ±q out to 
= ±^y^2m{h - A2/4eF)[l73. As expected (from 
the fact that these are different regimes of the same LO 
phase) the above outlined spectrum is in a qualitative 
agreement with that near hc2 , found in the previous sub- 
section. 



C. 



Fermion-Goldstone modes coupling in the 
Larkin-Ovchinnikov state 




FIG. 14: Lowest order Feynman diagram of a fermionic con- 
tribution, that corrects the superfluid stiffness. 



FIG. 15: Lowest order Feynman diagram of a fermionic con- 
tribution, that corrects the smectic compressional modulus. 



As discussed above, the LO state is characterized by a 
simultaneous presence of gapless Goldstone modes (par- 
ticularly "soft" for the fully rotationally invariant geom- 
etry of an isotropic trap) and gapless fermions. Conse- 
quently, for a complete description a coupling between 
these must be included. These can in principle be de- 
rived from the microscopic BCS Hamiltonian evaluated 
inside the LO state beyond a conventional mean-field 
treatment. Alternatively, the fermion- Goldstone mode 
couplings can be simply deduced based on symmetry con- 
siderations. The leading ones include, the coupling of the 
supercurrent and of the LO phonon to the quasi-particle 
current and number density, and are given by 



H. 

Ha 



P 



V0 • V^iV^ -I- h.c. (204a) 
(V0)2VV (204b) 
■ip^il; + (Vu • '^hVil)f + h.c. 

(204c) 



Because atom number conservation is effectively "bro- 
ken" in the superfluid LO state, above couplings must 
be supplemented by the anomalous (number violating) 
operators, such as e.g., {\7 (j))'^ tptp + h.c. Effects of these 
interactions on the fermionic and collective bosonic (0, u) 
spectral functions requires a detailed analysis, that par- 

180l-[l85l . 



allels studies of gauge fi elds 
[ll7j | and critical modes ll86 



Goldstone modes 

187l | coupling to electrons 
in gapless superconductors and metals. A preliminary 
analysis suggests that in the LO state these derivative 
couplings (enforced by the underlying gauge and spa- 
tial symmetries) only lead to a finite renormalization of 
model's parameters as well as Landau- like damping of the 
Goldstone modes. In addition to the above symmetry- 
dictated couplings, the prese nce of gapless fermions can 
generate Berry's phase terms [l55j |. that can qualitatively 
modify the conventional LO phonon dynamics derived in 
Sec lIIIDl We leave a detailed study of these interesting 
questions to future research. 



IX. LARKIN OVCHINNIKOV STATES IN A 
TRAP 



FIG. 16: Lowest order subdominant Landau-damping contri- 
bution to the superfluid stiffness coming from finite fermion 
density. 



The primary experimental application of our results 
is to polarized paired superfiuidity in trapped degener- 
ate atomic gases. It is thus crucial to extend our bulk 
analysis to take into account the effect of the trapping 
potential Vt{r), that in a typical experiment is well- 
approximated by a harmonic-oscillator potential. While 
a full analysis of the effect of the trap is beyond the scope 
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of this manuscript, in the present section we study this 
problem within the well-known local density approxima- 
tion (LDA). We note that several recent studies (e.g., 
Refs. % lii, MM, M, IM!) have also addressed polar- 
ized superfluidity in a trap. 

A. The local density approximation 

Much like the WKB approximation, LDA corresponds 
to using expressions for the bulk system, but with an 
effective local chemical potential /j,(r) = /i — Vt{r) in 
place of /i. The validity of the LDA approximation re- 
lies on the smoothness of the trap potential, with the 
criterion that Vt{r) varies slowly on the scale of the 
longest physical length A (the Fermi wavelength, scat- 
tering length, effective range, etc.) in the problem, i.e., 
{X/Vt{r))dVt{r)/dr <^ 1. Its accuracy can be equiva- 
lently controlled by a small parameter that is the ratio of 
the single particle trap level spacing SE to the smallest 
characteristic energy Ec of the studied phenomenon (e.g, 
the chemical potential, condensation energy, etc.), by re- 
quiring SE/Ec <C 1. Within the LO state the longest 
length is clearly the LO period, a = 27r/(7o, that near hc2 
is bounded by the coherence length (that near unitarity 
can be as short as inter-atomic spacing ~ R/N^^^, where 
R is the trapped condensate radius and N is the total 
number of atoms), and thus <^ R. Thus, in this regime, 
away from hd (where for a continuous BCS-LO transi- 
tion the period a is expected to diverge) the effects of the 
trap can be safely treated within the LDA. 

The generalization of the model of an imbalanced res- 
onant Fermi gas to a trap is straightforward: 

H= /'d3,Y^tzY!^^ + (y^(r)_^,(r))^t^, 
J \ 2m 

+ g^i^li^^i^^), (205) 

where V'<T(r) is a fermionic field operator with a Fourier 
transform Cko-- Henceforth, to be concrete, we shall focus 
on an isotropic harmonic trap (although this simplifica- 
tion can easily be relaxed) with 

Vt{r) = imw,V2, (206a) 

^ A^^, (206b) 

latter expression defining the cloud size R. Within 
LDA (valid for a sufficiently smooth trap potential Vt{r), 
see above), locally the system is taken to be well- 
approximated as uniform, but with a local chemical po- 
tential given by 

fi{r) = fi~^muj^r^, (207a) 



where the constant fi is the true chemical potential (a 
Lagrange multiplier) enforcing the total atom number 
N. The spatially-varying spin-up and spin-down local 
chemical potentials are then: 

fi^{r) = fi{r) + h, (208a) 
filir) = fi{r) - h, (208b) 

with the chemical potential difference h uniform. 

Consequently, within LDA the system's energy den- 
sity is approximated by that of a uniform system, (jl7l) . 
with the spatial dependence (via the trap) entering only 
through fj.{r). The ground state energy is then simply a 
volume integral of this energy density. Thus, the phase 
behavior of a uniform system as a function of chemi- 
cal potential, /x, translates into a spatial cloud profile 
through /i(r), with the critical phase boundaries /ic corre- 
sponding to critical radii defined by fic = IJ-{rc, h). 0, 113] 
As first predicted [il], this leads to a shell- like cloud 
structure that has subsequently been observed experi- 
mentally [34 37]. 

Within the LDA, we can furthermore deduce the effects 
of the trap on the structure of the LO state. We turn to 
this analysis next. 




FIG. 17: The collinear Larkin-Ovchinnikov state, compressed 
by an isotropic trap, Vt(r) (on the right), showing the en- 
hanced imbalance (P ~ —dzu) confined to the edge of the 
trap. 




FIG. 18: Possible alternative forms of the LO state confined 
in a tight isotropic trap. The form on the right trades oH' the 
cost of the curvature energy for the dislocation energy. 
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B. Trap-induced elastic distortion 

As discussed earlier, the Larkin-Ovchinnikov state 
breaks translational synimetry and as other crystals 
thereby exhibits a rigidity to stress. Confinement in a 
trap induces a variation in the chemical potential that 
(as observed experimentally in the BCS superfluid [13- 
ISTj l) has the effect of expelling the imbalance (polariza- 
tion) to the edge of the cloud. Since in the LO state 
the BCS order parameter domain walls "carry" the im- 
balance (the imbalanced fermions are confined to domain 
walls in Alo(-2), where their cost is minimum), in the LO 
state the increased imbalance is accommodated by the in- 
crease in the density of domain walls. Thus, the trap con- 
finement introduces an effective longitudinal stress that 
acts to compress the LO state (increase its local wavevec- 
tor go(r)), where the chemical potential is reduced. This 
tendency is captured by introducing a local longitudinal 
stress, a{r) into the LO phonon Hamiltonian 



(209) 

Above trap-phonon coupling arises by repeating our 
derivation in Secs lII|IlTl with the trap potential Vt{r), 
treating it within the LDA. The effective local chemical 
potential /Lt(r) enters through many contributions, but 
the leading order effect comes from a position-dependent 
wavevector generalization of Eq. (j25p 



qoir) = 2a 



vpir) ' 

2fi 



go 1 + 



(210a) 
(210b) 
(210c) 



In above, consistent with the LDA applicability we uti- 
lized the Vt{r)/iJL <ti 1 limit. Since Sqo/qo « —dzU sa 
—a/B, we deduce that the local stress is given by 



cr(r) 



1 B 

2 Tio/i 



n{r)Vt{r), 



(211a) 
(211b) 



where ctq = B/2 is the trap induced stress scale, and 
we inserted an additional factor of n{r)/nQ to crudely 
account for the breakdown of LDA based analysis near 
the cloud edge, where atom density, n{r) vanishes at R. 
As expected on physical grounds, (T(r) is a function that 
vanishes at the trap center and peaked on its outer shell 
set by R. 

With this, the distortion uo(r) is given by a standard 
form 



u(r) = / dVG,(r-r')a(r'), 



(212) 



in terms of a Greens function Gz{'fc) = dzG{r) that is a 
derivative of the smectic Greens function, with 



G.(r) 



dPq iq^e^'^'^ 
(2^)3 Kq\ + Bql ' 

-1 



i7rB\\z\ 



(213a) 
(213b) 



From this the trap-induced imbalance distortion, SP{r) 
is straightforwardly computed 



SP{r) w ~Pdzu{r), 

= P J^ dzGz{r-r')a{r'), 



(214a) 
(214b) 



. ^^nirmir)^P^^l--^), (214c) 

with final unnecessary crude approximation by construc- 
tion of cr(r) and its coupling to u recovering the expected 
result. 



C. Trap-induced LO 3d-2d dimensional crossover 

A microscopically accurate account of the trap in- 
side the LO state, and in particular its coupling to the 
Goldstone modes with appropriate boundary conditions 
(beyond above kludge treatment) is not currently avail- 
able. Perhaps a numerical analysis (e.g., numerical so- 
lution of BdG equations or full quantum Monte Carlo 
simulation) (52l . \5&L [70| can provide the desired descrip- 
tion. 

However, lacking such first principles analysis, we are 
constrained to proceed phenomenologically, working di- 
rectly with the effective Goldstone mode theory of u. Ex- 
panding the LO state about the trap-distorted state uo{r) 
discussed above, the distortion du{r) is again governed 
by the smectic Hamiltonian (jl42b[) . Although, some of 
the effects of the trap will enter through the position de- 
pendent elastic moduli B(r),K{r), we expect that the 
leading effects of the trap are incorporated through the 
physically motivated boundary conditions on u(r). Be- 
cause atom density vanishes at the edge of the cloud, 
the LO inner shell is expected to be surrounded by an 
outer shell of a fully polarized normal cloud lacking any 
positional order. Thus, we expect Neumann boundary 
condition 



fit ■ Vm — 0, 



(215) 



{fit the unit normal to the cloud's boundary) to be the 
most appropriate one to supplement the Euler-Lagrange 
equation for u{r). However, the appearance of forth 
derivative along r , Eq. (jl42bp requires that this bound- 
ary condition be supplemented by additional ones on the 
domain-wall normal, V_lw. 

Given the cylindrical form of smectic elasticity, 
Eq. (jl42bp . analytically implementing such spherically 
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symmetric boundary conditions is quite challenging. 
However, because our goal here is more modest, a qual- 
itative understanding of the trap-induced dimensional 
crossover can be obtained by simply using cloud size scale 
R to cutoff long scales appearing in a bulk analysis. Tech- 
nically, our analysis amounts to instead working with pe- 
riodic boundary conditions and a cylindrical trap. The 
calculational convenience of such boundary conditions is 
that they do not modify the form of the eigenmodes (still 
Fourier modes), and enter only through a restriction on 
the allowed eigenvalues, — 2'Knz/Lz, q± = 2Trn±/L±, 
with nz^± G Z. Thus, length scales associated with the 
trap and finite cloud size, Lz,L±, crudely enter through 
the minimum allowed momentum eigenvalues (roughly 
set by the cloud size Rz.,±), and thereby capture the spa- 
tial extent of the lowest phonon eigenmode even for cor- 
rect boundary conditions. 

A key observation that emerges from such treatment is 
that because of the "infinite" anisotropy of the bulk smec- 
tic modes in Eq. (|142b[) . with z ^ r^/A (or equivalently 
Qz ~ Ag^), the dimensional crossover is qualitatively dif- 
ferent than in systems with more conventional, scaling- 
wise isotropic elasticity. Namely, examining smectic bulk 
propagator Gq = l/{Bq^ -\- Kq\), it is clear that un- 
less Lz ~ L\/\ ^ L±, the phonon fluctuations will 
not be controlled by bulk modes. That is, as illustrated 
in Fig ll9l for any reasonably shaped trap (even quite 



anisotropic one, other than an extremely anisotropic 
with Lz ~ L'^/X), fluctuations will be controlled by the 
d ~ 1 dimensional "zero" modes, uo(r±), that are uni- 
form along the z axis, representing compression-free LO 
undulations. Indeed this is allowed because of the ex- 
pected Neumann boundary condition on u(z,rj^), that 
allows a zero-energy cost rigid displacement along z of 
the LO domain- walls. 

To demonstrate above expectation, we account for 
these "zero" modes by using a generalized mode expan- 
sion of u(r) that explicitly includes the z-independent 
modes 



M(r) = 



1 



LzL\ 



J2 uiqz,qi_)e:^-'' + 



(216) 



In terms of these the Hamiltonian becomes 



H?, — 



^— J2 [iBq'z+Kqi)\u^\' + KSq^Mi\ulJ'], 

_L 

(217) 
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and can be used to calculate smectic phonon fluctua- 
tions. For example using H!^ and equipartition, local 
root-mean-squared (rms) fluctuations are given by 



E 



T 



■E 



T 



z^^ Bql + Kql LzLi 



KqV 



T 



T 



Bql + Kqi 



1 r T 
17z I 



T L 



■Ani^L^/a) 



Kql 
T Ll 



8n^K L 



^ for Lz < L\I\ 
I 



(218a) 
(218b) 
(218c) 
(218d) 



where in the second line we approximated mode sums 
by integrals in a standard way. Indeed as summarized 
by (|218d|) . by comparing the two (bulk and zero mode) 
contributions to Urmsi anticipated on general grounds 
above it is clear that for isotropic (and even highly 
anisotropic) traps, the "zero" mode second contribution 
dominates as long as Lz < L\/X. 



harmonic level is given by 



S,, 



drd 



id-l 



-{drUo) 
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(vi«o) 



(219) 



Applying this to a spherically symmetric trap with 
Lz ^ L± fa R <^ R'^/X, we conclude that indeed fluctua- 
tions are controlled by the "zero" -mode LO phonons. We 
thus conclude that all the correlation-function properties 
of the LO state will be even more anomalous, character- 
ized by a d — 1 dimensional zero-mode action, that at a 



More detailed implication of these observations, particu- 
larly in contexts of specific experimental geometries re- 
main to be explored. We hasten to add, however, that 
by equipartition, the finite temperature the LO thermo- 
dynamics will nevertheless be dominated by bulk modes, 
simply due to their bulk-to-surface dominance. 
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5. 




LO SF in isotropic trap 



* ^ ~ Agl < ^ 
"infinitely" anisotropic fluctuations 

FIG. 19: Illustration of the = G (z-independent) phonon 
modes, even in an isotropic trap (indicated by a blue circle) 
dominating over the (scaling- wise) anisotropic, ~ Aq^ bulk 
smectic modes (indicated by an ellipse). 



X. EXPERIMENTAL IMPLICATIONS 

There is a large number of experimentally observable 
effects, that emerge from our study. A comprehensive 
treatment of these requires further extensive studies, that 
lie outside the scope of the present manuscript. Here we 
simply sketch out a few of the most important experi- 
mental signatures of our predictions. 



A. Larkin-Ovchinnikov order parameter 

As discussed in the Introduction and in SecEl at finite 
temperature in the thermodynamic limit (see below for 
the discussion in the trap) isotropically-iia.'p'ped Larkin- 
Ovchinnikov phase is characterized by a vanishing aver- 
age LO order parameter, i.e., (Ar , o) = 0. Reminiscent 
of 2d superfluids and crystals [130, Il4l| . this is a reflec- 
tion of its enhanced thermal fluctuations. As with these 
well-known examples and other topological phases, this 
does not however imply that the state is unstable (at 
least not in 3d), but that it requires a flner character- 
ization (e.g., correlation functions, topological defects, 
etc.) beyond a simple Landau order parameter. One of 
the experimental implications is, that, consequently the 
leading nonzero Landau order parameter characterizing 
the LO state is the translationally- invariant "charge" -4 
(4-atom pairing) superconducting order parameter, Asc, 



introduced in Eq. (I70ap . Thus, in the presence of thermal 
fluctuations the LO phase corresponds to an exotic state 
in which the off-diagonal order is exhibited by pairs of 
Cooper pairs, i.e., a bound quartet of ato ms, rather than 
by the conventional 2-atom Cooper pairs jll9j |. In 2d and 
3d this higher order pairing is driven by arbitrarily low-T 
fluctuation, rather than by a flne-tuned attractive inter- 
action between Cooper pairs, and therefore akin to 2d 
superfluids and crystals has no mean-field description. 
While a direct experimental probe of A^c may be chal- 
lenging, enhanced fluctuations in the LO state can be 
directly observed through correlation functions to which 
we turn next. 



B. Momentum distribution function 

Probably the most striking signature of the LO state 
is the novel form of the Cooper-pair momentum distri- 
bution function, rik — (Aj^Ak). As for a conventional 
resonantly-paired superfluid, rik should be accessible by 
a detuning sweep (controlled by a magnetic field) that 
projects the finite-momentum LO Cooper pairs Ak on 
the BCS side onto tightly bound molecules on the BEC 
side of the Feshbach resonance j], and then observed 
through a standard time-of fiight imaging of the result- 
ing molecular condensate. In contrast to a conventional 
bosonic and BCS condensates (that in a trap display a 
single peak, associated with a condensation into a low- 
est trap state, or its interaction-swelled equivalent), the 
LO condensate is expected to display a spontaneous re- 
ciprocal lattice of condensate peaks associated with its 
periodic structure, Figl^ 

Whi le in mean- field approximation rik is predicted 0, 
[ssl . [Tl3 to re semble a superfluid in an optical periodic 
potential [ll, Il88j | , the spontaneous nature of its transla- 
tional and orientational symmetry breaking, in the pres- 
ence of thermal fluctuations lead to important qualitative 
distinctions. Although the full anisotropic form is quite 
complex and best evaluated numerically, the asymptotic 



form of rik*^ can be readily obtained analytically. As 
calculated in Sec|Vl n^'~' exhibit a power-law (algebraic) 
peaks around harmonics q„ of the ordering wavevector 
qo, replacing the mean-fleld (5- function Bragg peaks of 
bosons in a periodic potentialfll Il88l |: 



nqo\ 



for d = 3, (220) 



where for simplicity we specialized to k = fc^z, the 
form factor Uq^ is given after Eq. (|124p and rj = 
qQT/{8n^/BK). This form is reminiscent o f (l+l)d 
Luttinger liquids and two-dimensional crvstals |l36l4l38l 
Il4l| . and is a reflection of the quasi-long-range order of 
the nonzero-temperature 3d LO state. While physically 
quite distinct, n^'~' is mathematically closely related t o 
the structure function of a conventional smectic [l29lll34 | . 
Another feature of n^'~' is the absence of the k = 
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FIG. 20: The finite momentum pairing at go and divergent 
3d smectic phonon fluctuations in the LO state are predicted 
to be reflected in the Cooper-pair center-of-mass momentum 
distribution function, nk, displaying power-law Bragg peaks, 
characteristic of the spatial quasi-long-range order. We ex- 
pect these to be observable through the time-of-flight mea- 
surements. 



condensate peak, that qualitatively distinguishes the LO 
state from a suDersohd [87l - [90| . where crystalline and su- 
perfluid orders merely independently coexist. We expect 
these novel features to be the smoking gun for the LO 
state, in principle observable in time-of-flight imaging. 



C. Structure function 

The density-density correlations, conventionally mea- 
sured in the reciprocal space using x-ray or neutron scat- 
tering is another important quantity that can be exper- 
imentally probed. The simplest is the static structure 
function, that, using p(r) w |Aio(r)p is straightfor- 
wardly computed (see SeclV|: 



E 



A 



2g„ 



(221a) 
(221b) 



where for simplicity we evaluated it at q = q^z and 
^2q„ is a form factor. S^'-'{q) contrasts with n^*^ by 
its insensitivity to the off-diagonal (i.e., superfluid) or- 
der. In three dimensions it also displays quasi-Bragg 
peaks, but at twice the reciprocal lattice vectors, 2n(7o, 
with the An^rj fluctuation exponent, and just like conven- 
tional smectics[129, 130, 135] (unlike ?^^*^) does exhibit 
the qn=o — peak. 

After a projection onto the molecular BEG state, a 
density profile p(r) for a large atomic cloud should be 
measurable in situ. Computing a Fourier transform of 
its correlations then allows one to test S^'^ , above. This 



FIG. 21: The structure function, S(q) for the 3d LO state, 
displaying power-law (as opposed to (5-function) Bragg peaks, 
characteristic of the LO superfluid's spatial quasi-long-range 
order. 



way, in principle the dynamic structure function Sicw (a 
Fourier transform of two-time density correlation func- 
tion) should also be accessible. 

Alternatively, the dynamic structure functi on can b e 
directly measured using Bragg spectroscopy 157l4l59l | , 
with its static limit obtained by integrating over frequen- 
cies. Based on successful measure men t s of the Bogoli- 
ubov mode in bosonic condensates [HI, Il59l |. we expect 
that the peaks in the 5*^2 (associated with its poles) can 
be used to give the dispersion of the collective modes, the 
superfluid phase and LO phonon u. 



^0(k) = Xo'Vp^fci+plS, 



w«(k) 



Kk\ 
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(222a) 
(222b) 



with the corresponding lifetimes remaining to be deter- 
mined. 

Related to the dynamic structure function is the gener- 
alized density and polarization (spin-imbalance) response 
function to a periodic potential. The periodic nature of 
the LO state makes it possible to use a perturbation at 
half the wavelength of the LO state to resonantly ex- 
cite long-scale collective modes, such as the polarization- 
dipole mode, that can be readily detected in a trapped 
gas. This was pro posed and carefully investigated for a 
Id LO state in Ref ll89l and should also apply in the two 
and three dimensional cases considered here. We leave 
the analysis of these and related experiments to a future 
more detailed study. 



D. Trap effects 

As is clear from the discussion in Sec lIXI a detailed 
study of the effects of the trap is essential for an explicit 
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contact with experiments. By introducing system bound- 
aries, a trap modifies the structure of and fluctuations in 
the LO state, that, aw ay fr om hd (where LO period can 
potentially diverge [96i Il02j [) can be treated within LDA. 



1. Fluctuations 

The detailed phenomenology of the LO state in a trap 
depends on the nature of the boundary conditions, that 
still remain to be understood. However, for large atomic 
clouds with small surface to volume ratio, the boundary 
effects are expected to be a weak perturbation to the 
bulk. As discussed in Sec lIXi in this case its effects can 
be taken into account phenomenologically i n the spiri t of 
the finite size scaling of critical phenomena [l3Ct Il64| by 
introducing a long scale cutoff (set by the cloud size R) 
into our bulk quantities. For observables related to the 
soft LO phonon u, that is strongly affected by fiuctua- 
tions this introduces an experimentally observable sensi- 
tivity to the cloud size, R and by extension to the atom 
number TV, trapping frequency LUt and the explicit trap 

anisotropy cj^ / ivf . 

Referring to the analysis of Sees. IVIIXl that give 
phonon rms fluctuations for a gas in a "box" of size 
L± X Lz, we conclude that in a isotropic trap with cloud 
size R ~ L± = Lz, the phonon fluctuations are large but 
finite, in 3d given by 



InR/a, (223) 



where we approximated the phonons by Gaussian fiuc- 
tuations, neglecting weak nonlinear elastic effects of 
Sec lVB 2] (that can become important on scales R > 
C±z7 (|128|) ). and assumed Neumann boundary condi- 
tions that allow for the = "zero" modes. The 
crossover from the bulk to the "zero" modes domi- 
nated regimes takes place for a cloud size R > R^, — 

One consequence of this result is that in a trap the 
thermally averaged LO order parameter, (Aio(r))"'" = 
2Aq„ (i?) cos (qo • r) no longer vanishes (as it does in an 
infinite bulk system), but strongly depends on the cloud 
size. For Neumann boundary conditions, its thermally 
suppressed amplitude in 3d is given by 



'^T/iio^/^h^-R/Ra (^^y^^ ^ (224) 



with Rq sa QqT/K the scale beyond which the "zero" 
modes become particularly important. Thus through 
R{N,LOt) the order parameter amplitude is also expo- 
nentially sensitive to the atom number N and trap fre- 
quency ujt- We note that the exponential dependence 
on i? is a consequence of the dominance of the Qz — 
"zero" -modes (under Neumann boundary conditions 
on u) over the anisotropic smectic bulk modes in an 
isotropic trap. In contrast, under Dirichlet boundary 



conditions that exclude these "zero" -modes, the expo- 
nential factor is suppressed, and the LO amplitude is 
a weaker power-law function of the cloud size, remi- 
niscent of 2d xy-systems, such as a superfiuid confined 
to two-dimensions. We note, that Kosterlitz-Thouless 
phase fiuctuation p hysics ha s been reported in 2d trapped 
atomic superfluids jlQOl flQlj . despite the finite trap size. 
We therefore expect our predictions for strong fluctua- 
tions effects in the LO phase to also be experimentally 
accessible. 

Finally we note that given our prediction that gener- 
ically the LO state is sandwiched by its descendents 
phases {SF^, FLsm, FL^, FL* , etc, rather than a sim- 
ple vacuum) the analysis and implementation of bound- 
ary conditions for a trap is further significantly compli- 
cated, particularly near hd. 



2. Phase diagram through shell structure 

As discussed and analyzed in Sec lIXl an even stronger 
effect of the trap is that it leads to an effective lo- 
cally varying chemical potential ^{r) = /i(l — r"^ /R^), 
and therefore gives slices through the chemical potential 
phase diagram as a function of radius r. We thus ex- 
pect that the phase diagram in Fig|4] can be "imaged" in 
the spatial cloud profile, with critical phase boundaries 
/ic*'' (to phase i) translating to critical radii of shells de- 
fined by /ii*'' = fiircK h). Given the past success of such 
phase detection with bosons in optical potentials (ex- 
hibiting "wedding cake" profiles) and the uniformly 
paired superfluids and Fermi liquids of an imbalanced 
Fermi gas [5, '34*, 'ss', '3?] , we expect that similar identifl- 
cation will be possible for some of the LO liquid crystal 
phases. 



E. Response to rotation 

Consistent with its neutral superfluid order, above a 
critical rate of rotation, Vld a FFLO state responds to 
an imposed rotation by nucleating quantized vortices. 
Because of its crystalline superfluid form, even within a 
mean-fleld description its vort ices are p redicted to form 
a rich variety of vortex Iattices [l92l4l94j |. that depend on 
the nature of the "host" FFLO state. 

We predict a number of additional interesting vortex 
features that are special to strongly fluctuating LO state. 
One distinguishing feature of the smectic LO state stud- 
ied here is its uniaxial anisotropy, that is strongly tunable 
with species imbalance. It manifests itself in the (quanti- 
tative) anisotropy of the superfluid densities (stiffnesses) , 

pjr ^ pi, Eq. (15^ . and the qualitatively anisotropic dis- 
persion of the LO phonon, Eq. ((89)) . The former leads 
to vortices with a spontaneously elliptical vortex cores, 
pSSp . that we expect to form a centered rectangular 
lattice for vortices (rotation axis) oriented transversely 
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to qo- This striking feature should be easily identifi- 
able in the rotation experiments of the type previously 
used to identiiy phase separation in imbalanced Fermi 
gases 0, 113] • In this transverse geometry we also expect 
vortices to be pinned to LO phase fronts, localized in ze- 
ros of domain- walls in Alo(j), akin to vortices in layered 
superconductors e.g., BlSCCOflOS^. In contrast for ro- 
tation about an axis along qo, a conventional hexagonal 
vortex lattice is expected. 

Another striking feature of the LO state found in 
Sec I VII is the 7r-vortex bound to a half- a/2-dislocation, 
with a 1/4 energy cost (in thermodynamic limit) of a con- 
ventional 2tt vortex. We therefore predict that it is this 
fractional vortex that will be preferentially induced by 
the rotation. One experimental consequence of this is the 
reduction in the lower-critical frequency ilci down to a 
1/2 of the conventional value 2mR2 lni?/C in a large cloud 
of radius R. Concomitantly, in addition to its imbalance- 
tunable anisotropy (that diverges near the upper, hc2 
phase boundary), the vortex lattice is characterized by 
a vortex density = ^^f^j that is increased by a factor 
of 2 from the conventional value of a rotated paired super- 
fluid. Perhaps the most striking feature that we predict is 
that a rotation-induced lattice of half-dislocations must 
accompany and be locked to this 7r-vortex lattice. We 
expect above novel features to be most pronounced for a 
rotation axis transverse to qo, with a more conventional 
rotational response for $7 || qo. Finally, the LO state 
may exhibit more than one critical rotational velocity 
rici corresponding to distinct onsets of the penetration 
of half-integer and integer vortices. 



F. Fermionic excitations 

As we have seen above, in addition to the low- 
energy Goldstone modes, (f>,u, the LO state is charac- 
terized by gapless fermionic excitations associated with 
the imbalanced atom, that are localized on zeros of 
domain- walls in Aio(r). The most direct probe of 
these mid-gap fermionic states is through the Fcshbach 
resonance [196| and the RF spectroscopics! 149, 197 200], 
with the momentum- resolved extension pOll 12021 ] . allow- 
ing one to measure the dispersion and Fermi surfaces 
pockets of the Andreev states calculated in Sec lVIIll An 
observation of coexisting gapped and (Andreev) gapless 
features in these spectroscopies, along with the superfluid 
phase coherence in rik would provide strong evidence for 
realization of a FFLO state. 

Such measurements can be comp lemented with shot- 
noise correlation spectroscopy [20j|, that has been suc- 
cessfully used to probe bosonic Mott insulators flj and 
fermionic paired condensates [204] . In the LO state, we 
expect pairing shot-noise correlations to be peaked for 
-kp + qo/2,U — > + qo/2,t and - qo/2,i< — > 
— kp- — qo/2,t atom pairs. In contrast to a conventionally 
paired BCS state, shot-noise correlations are furthermore 
anisotropically distributed in the center-of-mass pair mo- 



mentum, with peaks around ±qo, reflecting the sponta- 
neous nematic anisotropy of the LO state. 

The low-energy excitations in the LO state can also be 
probed less directly through thermodynamics, modifying 
the power-law in T behavior of, for example, the heat- 
capacity. They will also manifest themselves in thermal 
transport, though its experimental implementation in the 
context of trapped atomic gases remains an open prob- 
lem. 

Probably the most direct way to detect the existence 
of the LO state is to simply image the population species 
imbalance. Because imbalanced fermions are confined to 
Andreev states localized on the LO domain-walls, the 
polarization density is periodic, with its phase locked to 
that of ALo(r). As illustrated for the thermally-averaged 
LO order parameter, in a trap the amplitude of this 
periodic component of polarization will be nonzero and 
strongly R dependent. 

Quite clearly, significant detailed theoretical analysis 
is necessary beyond above qualitative discussion, but is 
left for a future research. 



G. Phase transitions and novel phases 

In addition to our findings of the rich fluctuation- 
driven phenomenology of the LO (superfluid smectic) 
state, in Sec lVIII we predicted a variety of putative de- 
scendent phases, that emerge by disordering the LO state 
through a set of continuous phase transitions. If in- 
deed realized as stable phases (something that our phe- 
nomenological approach is unable to determine for any 
specific system) organized into phase diagram in FigUl 
then many of the predicted features should be readily 
detectable in experiments on imbalanced resonant Fermi 
gases. Each of these states {SF^ FLm, FLsm, etc.) ex- 
hibits its own qualitatively distinct phenomenology, dis- 
cussed in Sec lVIII where they were defined. For example, 
Bragg peaks in the time-of-flight images can distinguish 
the periodic SFsm (LO superfluid smectic) state from 
the homogeneous SFj^ (superfluid nematic), which are 
in turn distinguished from the FL^g^ and FLm (normal 
smectic and nematic) by their superfluid properties, bro- 
ken spatial symmetries (periodicity and anisotropy), col- 
lective modes, quantized vortices, and condensate peaks. 
Standard thermodynamic signatures (e.g., heat capacity) 
will identify the corresponding phase transitions, though 
some associated with strongly non-meanfield topological 
type between two gapped (disordered) phases (e.g., FL* 
and FL) maybe more difficult to detect. In a trap (see 
above) , the most vivid manifestation of these states is the 
appearance of a shell structure, corresponding to slices 
through the chemical grand-canonical phase diagram of 
an imbalanced resonant Fermi gas, Fig|31 Again, we leave 
the detailed analyses of all these features to extended fu- 
ture studies. 
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XI. OPEN QUESTIONS 

While this manuscript addresses a broad range of phe- 
nomena associated with the Larkin-Ovchinnikov state, it 
leaves many interesting questions open. Certainly the 
most important of these is the long-standing question of 
the range of energetic stability of the crystalline super- 
conductor discussed in the Introduction. If the state is 
indeed stable over a sufficiently broad range of detun- 
ing and imbalance to be experimentally accessible, is its 
lowest energy form indeed the striped collinear LO type, 
assumed throughout this manuscript? While for large 
atomic clouds and shallow traps (such that LDA remains 
valid) we expect only a small deformation of the LO state 
near the boundaries of the phase, for tighter traps a more 
detailed treatment of the trap is necessary, and may lead 
to a distinct global form of the LO state, such as the 
"onion" and "radial" structures illustrated in Fig.([T8|. 
To address such questions undoubtedly requires numeri- 
cal solutions in experiment-specific geometries. 

Furthermore, the nature of the (2d and 3d) tran- 
sition into the LO state at the lower-critical Zeeman 
field /ici,!!^, and the extent to which it resembles a 
commensurate-incommensurate transition (as in ld[95!, 
[l08]) remains an open question. More broadly, in Sec I VIII 
we predicted a number of novel LO descendent phases ad- 
jacent to the LO smectic superfluid state, but have left 
open the detailed nature of their phenomenology, sta- 
bility to quantum and thermal fluctuations, as well the 
nature of the associated phase transitions. Similarly to 
the LO state, these phases are expected to exhibit gapless 
fermionic excitations coupled to their Goldstone modes. 
Understanding the effects of these fermionic modes on the 
properties of the phases and the associated transitions re- 
main wide open and extremely interesting problems. 

Finally, as is clear from the discussion of the previous 
section, much detailed theoretical analysis remains to be 
done to make contact of our general predictions with spe- 
cific experiments. We leave these and many other inter- 
esting questions to future studies. 



XII. SUMMARY AND CONCLUSIONS 

To summarize, we studied a wide range of fiuctuation 
phenomena in a LO state, expected to be realizable in 
an imbalanced resonant Fermi gas. Starting with a mi- 
croscopic description of a resonant Fermi gas, supported 
by robust model-independent and very general symme- 
try arguments we have demonstrated that in an isotropic 
trap the LO state is a gapless superfluid smectic liquid 
crystal, whose elastic moduli and superfluid stiffness we 
derived near hc2- Consequently, the state is extremely 
sensitive to thermal fiuctuations that destroy its long- 
range positional order even in three dimensions, replacing 
it by a quasi-longer range order of the resulting algebraic 
quantum smectic state, characterized by power-law cor- 



relations akin to a system tuned to a critical point or 
two-dimensional xy-model systems. We showed that this 
exotic state also exhibits vortex fractionalization, where 
the basic superfluid vortex is half the strength of a vor- 
tex in a regular paired condensate, and is accompanied by 
half-dislocations in the LO smectic (layered) structure. 

Studying the fluctuation-driven disordering of the LO 
smectic, we predicted a rich variety of descendant phases 
such as the superfluid {SFjy) and Fermi liquid nemat- 
ics (FLn) and the fractionalized nonsuperfluid states 
(FL*), that generically intervene between the LO state 
and the conventional BCS superfluid (at low population 
imbalance) and a conventional Fermi liquid (at high pop- 
ulation imbalance). We outlined a large variety of ex- 
perimental implications of our flndings, but leave their 
detailed analysis to future studies. 
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Appendix A: Ginzburg-Landau expansion 

In this appendix we provide some of the technical 
details necessary to derive the quartic interaction in 
Aq of the Ginzburg-Landau expansion appearing in 
Eas. ([30| . ((50)) and in particular the current-current con- 
tribution (jin]), that determines the transverse superfluid 
stiffness p-^. As outlined in the main text, we use the 
coherent-state imaginary-time path-integral formulation 
of the BCS partition function. We decouple the quar- 
tic fermion interaction g by introducing the Cooper-pair 
Hubbard-Stratonovich field A(x), and formally integrate 
out the fermionic atoms, obtaining 



(Al) 



where the effective Ginzburg-Landau action is given by 



5e//[A*,A] = -In 



(A2) 



where Sr is the microscopic action defined in Eq. (j33l) . 

Taylor-expanding in powers of Sint, Eq. p5bl) we 
obtain 
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^e//[A,A] = -InZo-ln 



(A3a) 
(A3b) 



where 



-So 



e ^"/Zq. Above, all odd-power in A 



terms clearly vanish, and the quadratic A*eA term has 
already been analyzed in the main text [1, 13, [t^ ■ 

Thus, we focus on the contribution quartic in A, that 
is a connected forth cumulant. Using the definition of 
Sint, Eq. (I35b|) we find 



4! 



4 \ c 
■mi/0' 



(A4a) 
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F(xi,X2,X3,X4)A;^Ax,A;3A, 



(A4b) 



where 



with the noninteracting fermionic Green's function (in 
Fourier space) as usual given by 



1 



(A6a) 
(A6b) 



y(xi,X2,X3,X4) = (A5) 
G?(X2 - X1)GJJ(X2 - X3)G?(X4 - X3)G;(x4 - Xi), 



Because we are interested in r-independent A(r), the t 
integrals can be taken, giving Si = J drHi, with 
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Hi = -j (2^)'^(5'*(qi-q2+q3-q4) 



^ y(qi +qi,qi +q2,qi +q3,qi +q4)A*^(qi)Aq,(q2)A*^(q3)Aqi(q4) 
qi=±q 



-4l^(q + qi, q + q2, -q + qs, -q + q4)A;(qi)Aq(q2)Alq(q3)A_q(q4) 



(A7) 



A Taylor-expansion of F(q„i + qi,q„2 + q2,q«3 + 
qs, qn4 + q4) in q^ gives Hi = Hf^ + where Hf'' 

has already been computed by LO in their mean-field 
approximation ^TTj and in Ref. from the expansion of 
the BdG ground state energy of the FF state discussed 
in SecHm 

We focus on the second H\ term, and Taylor expand 
it to second order in q^ 



F(q + qi , q + q2, -q + qs, -q + q4 

{quqij + quqzj + qiiqi] + q^mj) 



(0) 



(A8a) 



(A8b) 



that gives the key current-current interaction vertex, 
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(2) 2 



-I 



rlf.irl'^Jc ~ ^ 

^^G°{k, a;)2a,G°(q - k, -a;)a,G°(-q - k, -uj), 
(q-k),(q + k)j- 



(q-fcFk)i(q + fcFk)j 



(k+q)^ 



u},£,k {icj - s + hy{iijj + E- vpk ■ + ^ + + £ + vp^ " 1 + ^ + ^) 



We first carry out the uj integral using simple identity, 
duj 1 



2n {iuj — ay{iuj — 6)^(ia; — c)^ 



g3 



cku 



1 



dadbdc J 27r {iw — a){iuj — b){iLu — c) 



Q3 



dadbdc 



e(-a) 



+ 



e{-b) 



+ 



e(-c) 



{a — b){a — c) {b — c){b — a) (c— a)(c — 6) 



S{a) S{b) 6{c) 
b'^c^ c^a? a^b"^ 

N 2a — b — c ,x 2b — c — a 

-2B(-«)7:-^:t^7:— To - 2e(-6)- 



'(a-6)3(a-c)3 
wliicli, when used in the expression above, gives 

4{2s + £^)e{h-e) 



(6-c)3(6-a)3 



- 2e(-c) 



2c - a - 6 
(c-a)3(c-6)3' 



(A9a) 
(A9b) 

(A9c) 



(AlO) 



(fcij-k- q)j(fcFk + q)j 



e,k 



+ 



(2£+^-3z;Fk.q)e(e+^ 



2 ^ 

ft. — vpk ■ q) 



+ ^ + ■ q)Q(g + ^ + + ^'f k ■ q) 
4(t.Fk-q)3(2£+£+z;^^k-q)3 

(5(e- ft) 



4(t>i.k-q)3(2£- 



r_ 

2m 



Vpii ■ q)3 



[(2ft- 



2m' 



(t'Fk-q)^]- 



+ - 



+ + • q) 5{e+ ■§- + h-VF\^- q) 



+ 



AivF^k ■ q)2(2ft + £ + i;j.k • q)2 4(i;f k • q)2(2ft + £ - t;^k ■ q)2 
(fcFk-q)j(fcFk + q)j 



[(2/i+£)2-(t;Fk.q)2 



(All) 



In above we split the k integration into an integra- 
tion over its orientations k and magnitude k, and ap- 
proximated the latter using a nearly constant density of 
states at the Fermi energy, equivalently, ignoring small 
quadratic contributions in ^fc = fc — fc^, valid in the BCS 
limit. However, in contrast to a standard isotropic cal- 
culation, here integral over k is nontrivial because of the 
anisotropy introduced by q. That is, we used 



where the density of states per spin is 

,3/2 



N{e) 



m 



.1/2 



2l/27r2ft3^ 

Nie + ep) 
AT/ \ 3 n 



ce 



1/2 



f 2) 

By symmetry w-^- is clearly uniaxial along q. Thus it 
can be written as 



(A13a) 



(A13b) 



dkk^^jdn,..., 

deN{e)^ Jdn^..., 
Niep) J de^ j dVl^..., (A12) 



where 



(2) 



51 

92 



1 (2) 1 (2), , 

2< - 2^^^^^'*^' 



1 



.(2) 



O (2). , 



(A14) 

(A15a) 

(A15b) 
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with the sums over repeated indices implied. Using this 
inside Eq. (|A14p we find: 



51 = lN{ep)P^l- 



l-(k.q)2 



_ iV(ep)fc|, 
N{eF)k 



da — 

[a2 



1-^2 



24<z4 



-ai 



2/i 



(Al6a) 
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3(k . q) = 



1 



iV(eF)fc| 
24^4 

^ jV(ef)fc| 

,4 „4 



da — 

\a' 



2h 



[(2/i4 

3(7^ - 1 



2m' 



4fe^ 1 2 



„ 4 4 -a2 I I , (A16b) 

24% \VFqoJ 

where in above we evaluated q at the dominant dispersion 
minimum value of Qq and ignored subdominant q^/kp 
terms. 

The dimensionless functions ai{x) and a2{x) defined 
by the above polar angle (a) integrals are given by 



ai(a;) = 







da- 
1 

2^ 



1 



[a2 - a;2] ' 



1 + 



In- 



1 



4x3 -1-.' (^1^^) 
2.07, at hc2 {xc2 = l/ac2 = 5/6) transition, 



a2{x) 



da- 



3a' 



1 



[a^ - x^] 



X — 3x 



x^ 
i(3a;4 



2x' 



l)lni 



l+x 



2x3(1 -.2) - (A17b) 

~ —5.75, at hc2 {xc2 1/q?c2 = 5/6) transition, 

with the final values obtained by their evaluation at the 
upper-critical field hc2- Although naively the expressions 
appear to be linearly divergent, a careful analysis of the 
integrals regularized with an infinitesimal imaginary part 
(associated with the analytic continuation to real fre- 
quencies, iu! ^ uj + jO+) give well-defined values given 
above. 



Appendix B: Fluctuation in a finite "box" 



In the main body of the paper, Sec lIX Cl we demon- 
strated that due to the smectic phonon dispersion 
anisotropy, for Neumann boundary conditions the zero 
qz = mode dominates fluctuations over the bulk Qz ~ 
Xq^ modes and must be explicitly taken into account. 
This leads to an effective dimensional reduction for a sys- 
tem confined in any reasonably isotropic trap. Here, we 
rederive this result within a complementary, fully contin- 
uum description and contrast the behavior in the smectic 
state with that of the more familiar isotropic xy-model. 



1. Isotropic models 

For an xy-model (or really any model with isotropic 
spatial dispersion) described by 

H=IJ dU^u)^, (Bl) 

the fluctuations inside a "box" of aspect ratio Lz x Lj_, 
with Lz < i_L, for free (Neumann) boundary conditions 
are given by 



(«'(x)> = 



dqzd^-^qi^ 1 



1 



{2nr ql + qV 



(B2a) 



—L^^'^ > L'i-", for d<2,Lz<L^, (B2b) 
Liz. 



where the first term is due to the qz — ^ zero-mode and 
is larger than the second bulk modes contribution when 
< -^^_L- It accounts for fluctuations in the reduced 
dimensionality of a d — 1 dimensional " film" of thickness 
Lz ■ We note that for an isotropic dispersion the crossover 
to a lower-dimensional scaling takes place only when the 
actual geometrical aspect ratio of the system is "film"- 
like, i.e., anisotropic. 



2. Smectic models 

We contrast the above standard finite-size scaling with 
that of a model with a smectic dispersion, described by 



2 



d'^X {{dzuf + {y\uf 



(B3) 



In a finite "box" , its fluctuations are instead given by 



(u2(x)) = 



dqzd 



<d-i 



{2nr ql + q\' 
^L^-^ + Lf-^y^ 



(B4a) 



^ l^^i"' » Lf-''^'^ for d < 3, < Ll. 
Lz 

(B4b) 

Now, clearly the first zero-mode (g^ = 0) term domi- 
nates over the bulk contribution (second term) for Lz < 
L\ ^ L_i_, and does so even in an isotropic box with all 
dimensions L. Thus, a smectic confined in a geomet- 
rically isotropic environment is effectively deep in the 
lower, {d — l)-dimensional "film" regime with phonon 



fluctuations scaling like 



» ^(3-d)/2^ 
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